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ABSTRACT 

We present 4 square degree weak gravitational lensing survey of subhalos in the nearby Coma cluster 
using Subaru/Suprime-Cam. Thanks to large apparent size of cluster subhalos, we detected in a model- 
independent way 32 subhalos down to the order of lO^'^ of the virial mass of the cluster. Weak-lensing 
mass measurement of shear-selected subhalos allows us to investigate subhalo properties and the 
correlation between subhalo masses and galaxy luminosities, for the first time. The mean distortion 
profiles stacked over subhalos show a sharply truncated feature which is well-fitted by Navarro Frcnk 
& White (NEW) mass model with the truncation radius, as expected by the tidal destruction by the 
main cluster. We also found that subhalo masses, truncation radii and mass-to-light ratios decrease 
toward the cluster center. The subhalo mass function, dn/dlnMsuh-, in the range of two orders of 
magnitude in mass, is well described by a single power law or a Schcchter function. Best-fit power 
indices of l.llj^Q'jg for the former model and 0.94J^Q2g for the latter one are in remarkable agreements 
with slopes of ^ 0.9 — 1.0 predicted by the cold dark matter paradigm. The tangential distortion 
signals in the radial range of 0.02 — 2 h~^Mpc from the cluster center show a complex structure which 
is well described by a composition of three mass components of subhalos, NEW mass distribution as 
a smooth component of the main cluster, and lensing model from large scale structure behind the 
cluster. Although the lensing signals are one order lower than those for clusters at z ^ 0.2, the total 
signal-to-noisc ratio, S/N — 13.3, is comparable to or higher because a large number of background 
source galaxies compensate low lensing efficiency of the nearby cluster. 

Subject headings: galaxies: clusters: individual: Coma Cluster (A1656) - gravitational lensing: weak 
- X-rays: galaxies: clusters 



1. INTRODUCTION 

The cold dark matter (CDM) concordance cosmology 
has had considerable success in explaining various obser- 
vations on large scale, such as cosmic microwave back- 
ground (Komatsu et al. 2011; Hinshaw et al. 2012). It 
gives us the initial condition for the hierarchical struc- 
ture formation, involving in mass assembly histories of 
halos, for high-resolution A^-body simulations and ana- 
lytical models. In hierarchical clustering, less massive 
halos are accreted into more massive halos, subsequently 
eroded by effects combined with tidal stripping and dy- 
namical friction of the host halo, and eventually be- 
come its smooth component. Since galaxy clusters are 
the most massive virialized object in the Universe, the 
central regions of subhalos will survive under the over- 
density field until the recent epoch and constitute their 
population. The CDM numerical simulations and ana- 
lytical approaches predict that subhalo mass functions 
at the intermediate and low mass scales follows a power 
law, dn/dlnMsnh oc M^"^ with slopes of ~ 0.9— 1.0 (e.g. 
Taylor & Babul 2004, 2005a,b; Oguri & Lee 2004; van 
den Bosch et al. 2005; Diemand et al. 2004; De Lucia 
et al. 2004; Gao et al. 2004b; Shaw et al. 2006; Angulo 
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et al. 2009; Giocoh et al. 2010; Klypin et al. 2011; Gao 
et al. 2012). 

Observations of cluster subhalo properties, such as 
mass function and spatial distribution, provide us with 
a deeper understanding of its mass assembly history and 
make a stringiest test of CDM predictions on scales of 
less than several Mpc. A characteristic feature of sub- 
halo mass function is also critically important to con- 
strain the nature of dark matter, because it depends on 
the particle mass of dark matter. Furthermore, a study 
of the correlation between galaxy properties and subhalo 
masses, incorporating different data-sets, shed important 
insights on the physics of galaxy evolution associated 
with dark matter. Thus it is of paramount importance to 
measure the mass function directly from the observation 
without any assumptions on the relation between dark 
matter and luminous matter and dynamical state of the 
system. It is difficult, though, to infer subhalos' masses 
from visible matter, such as galaxies, because relevant as- 
sumptions on mass distribution extending beyond galax- 
ies and dynamical state of the objects are required. Here 
a weak gravitational lensing analysis plays an important 
role. Weak lensing analysis measures a coherent dis- 
tortion pattern of background galaxy images caused by 
gravitational field of the system and thus avoids any as- 
sumptions mentioned above (e.g. Bartclmann & Schnei- 
der 2001). However weak-lensing signal is obtained by 
averaging over a certain number of background galaxies 
and thus only the mass information over a scale of several 
arcminutes is obtained. Previous weak lensing studies or 
joint strong- and weak- lensing studies (Broadhurst et al. 
2005; Okabe et al. 2010b, 2011, 2013; Oguri et al. 2010, 
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2012; Umctsu & Broadhurst 2008; Umetsu et al. 2010, 
2011; Applegate et al. 2012; Hoekstra et al. 2012) mainly 
focused on clusters in the rcdshift higher than 0.15 or 
so, because of good lensing efficiency and high quality 
imaging data taken by the wide-field camera mounted on 
ground-based telescopes, such as the Subaru Prime Fo- 
cus Camera (Suprime-Cam; Miyazaki et al. 2002) on the 
8.2m Subaru Telescope. However weak lensing analysis of 
those cluster at z ~ 0.2 is very difficult to detect subhalos 
of which masses are at the order of 10^^/i~^Mq, because 
the apparent resolution ( ^ 0.'2 — 0.'5 or less) is too low 
to detect the truncation feature in lensing signal. Stack- 
ing lensing studies for member galaxies help to improve 
this disadvantage (Natarajan & Springel 2004; Natarajan 
et al. 2007). The staking method computes lensing dis- 
tortion signals centering around member galaxies, and 
thus is independent of subhalo selections and increases 
the signal-to-noise ratio in lensing signals thanks to large 
number of sample. Hence, the stacking analysis enables 
us to measure a mean mass of subhalos corresponding to 
member galaxies. However, when we derive subhalo mass 
function and/or conduct statistical study of subhalos, we 
cannot avoid the assumption between the luminosity and 
mass for subhalos. 

Weak-lensing study of nearby clusters (z ^ 0.1) over- 
comes all these problems mentioned above. Contrary to 
weak-lensing studies of clusters at z ~ 0.2, there are three 
great advantages for the analysis of dark matter subha- 
los: First, a large apparent size enables us to easily re- 
solve less massive subhalos inside clusters. Second, sub- 
halos are separated from main cluster center and other 
subhalos enough to ignore their lensing contamination 
in subhalo mass measurements. Third, a large angular 
scale provides a correspondingly large number of back- 
ground galaxies, which leads to low statistical errors and 
compensates low lensing efficiency to achieve high signal- 
to-noise ratio. The last one also plays an important role 
in cluster mass measurement. For example. Coma clus- 
ter is at redshift Zc = 0.0236, giving a large apparent 
size which is ~ 7 times larger than those of clusters at 
z ~ 0.2. We thus use the area of 50— 100 square minutes 
or more in weak-lensing mass measurements of subha- 
los of which masses are larger than ^ 10^^h~^MQ. In- 
deed, Okabe et al. (2010a) has demonstrated the power 
of weak-lensing analysis of Coma cluster and discovered 
less massive subhalos. 

Here we report the result of 4.1 square degree weak 
gravitational lensing survey of subhalos in Coma clus- 
ter by 18 pointing observations {Re and V bands) 
with Subaru/Suprime-Cam to directly measure subha- 
los' masses and their mass function. This paper is a 
continuation of our previous work (Okabe et al. 2010a). 
Okabe et al. (2010a) used archival Subaru/Suprime-Cam 
data with Re band. The imaging data cover the central 
and the south-west regions (two pointings; see also Fig- 
ure 1), of which area is in total ~ 0.5 square degree. On 
the other hand, our new data significantly improve the 
quality for weak-lensing analysis. First, the data covers 
out to the cluster outskirts, which enables us to study 
the radial dependence of subhalos properties. Area frac- 
tions within for previous and new data account respec- 
tively for ^ 10% and ^ 80% within r2oo inside of which 
the mean interior density is 200 times the critical mass 
density at cluster redshift. Second, the exposure time is 



deeper than the effective one for previous weak-lensing 
analysis (~ 16 min;Table 1), which increases the number 
of background galaxies and thus suppresses the noise on 
intrinsic ellipticity. Third, we use two filters to make a 
secure background galaxies to avoid a contamination of 
unlensed member/foreground galaxies in our shear cat- 
alog. Therefore, new data enables us to conduct a sys- 
tematic survey for cluster subhalos, for the first time. 
We describe the details of data analysis in Sec 2, in- 
cluding shape measurements, background and member 
selections and modeling of background lensing signals. 
In Sec. 3, we define subhalos from lensing signals, mea- 
sure model-independent projected mass, conduct stacked 
lensing analyses and evaluate systematic errors including 
the purity of the sample. We also present a galaxy-galaxy 
lensing study for luminous member galaxies in Sec. 4, 
which is complementary to and independent of the anal- 
yses in Sec. 3. We describe the main cluster mass mea- 
surement in Sec. 5. Finally, we discuss subhalo mass 
function, subhalo properties, and future studies in Sec 6 
and state our conclusion in Sec 7, respectively. Through- 
out this paper, we use the cosmology of rimo = 0.27, 
fl\ = 0.73 and Hq = lOOh km s^^ Mpc~^. One arcmin 
corresponds to 20 h~^kpc. 

2. DATA ANALYSIS 

2.1. Survey Observation and Image Processing 

We observe Coma cluster with Suprime-cam (Miyazaki 
et al. 2002) at the Subaru 8.2-m telescope, in Re and V 
bands, through March and April, 2011. The Re band 
data is used for the wide-field weak lensing analysis, and, 
a combination with V band data is used to minimize 
a contamination of member and foreground galaxies in 
shear catalog. The survey is covered by a mosaic of 18 
pointings, namely, comalO ... coma43, as shown in Fig- 
ure 1. The overlapping length between each pointing is 
2 arcmin. The survey area is in total 4.1 sqdeg. A maxi- 
mum projected radius from the brightest cluster galaxy, 
NGC 4874, reach ^ lOOf which is comparable to the 
cluster virial radius rvir (Sec. 5). The Re band data 
of comaSO was retaken due to the low number of back- 
ground galaxies. The typical exposure times for Re and 
V bands are 24.5 and ^ 14.0 min (Table 1). We also use 
two Re imaging data retrieved from the Subaru archival 
data (SMOKA^). 

We use the standard pipeline reduction software for 
Suprime-Cam, SDFRED (Yagi et al. 2002; Ouchi et al. 
2004) modified for new CCD, for flat-fielding, instru- 
mental distortion correction, differential refraction, PSF 
matching, sky subtraction and stacking. The seeing for 
each pointing is shown in Table 1. An astrometric cali- 
bration was performed using point sources of the 2MASS 
catalog (Skrutskie et al. 2006). The typical residual val- 
ues is no larger than the CCD pixel size. Photometric 
calibration was carried out by fitting point sources de- 
tected in each data with stars of SDSS DR8 photometry 
(Eisenstein et al. 2011), taking into account the differ- 
ence between their sensitivities. The archival data with 
old CCDs is reduced by the same procedure using the 
SDFRED for old CCDs. 

2.2. Weak Lensing Distortion Analysis 
^ http://smoka.nao.ac.jp/index.jsp 
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Fig. 1. — Cadence design for Coma cluster subhalo survey. The 
horizontal and vertical axes are R.A. and Dec with an offset dis- 
tance from NGC 4874, in unit of arcmin, respectively. The data 
name (Table 1) is shown in the middle of each pointing. The over- 
lapping length between pointings is 2 arcmin. The survey area is 
in total 4.1 sqdeg. Hatched regions are the areas of archival data 
(core and sub). The summary of imaging data is described in Table 
1. 



TABLE 1 

Imaging data 



Our weak Icnsing measurement follows Kaiser et al. 
(1995), so-called KSB-H method, using the IMCAT pack- 
age ^ with some modifications in a similar way as Umetsu 
et al. (2010) and Oguri et al. (2012). We measure the im- 
age ellipticity from the weighted quadrupole moments of 
the surface brightness of objects detected in each Re band 
imaging data (Table 1). The PSF anisotropic correction 
is conducted with the same way as Okabe et al. (2010a, b, 
2011). We select bright unsaturated stars in the half- 
light radius, r/i, and magnitude plane to estimate the 

stellar anisotropy kernel, g* = (J^s*m)a^s*' where P^ 
is the smear polarisability matrix, Ca is the image ellip- 
ticity, and the quantities with asterisk denote those for 
stellar objects. Following KSB method, PSF anisotropy 
is corrected with the equation 



pa/3 * 



(1) 



We estimate q"{9) at each galaxy position, 0, using a fit- 
ting function of a second-order bi-polynomials of vector 6 
with iterative cr-clipping rejection. The data region is di- 
vided into several rectangular blocks based on the typical 
coherent scale of the measured PSF anisotropy pattern. 
The median for two components of residual stellar ellip- 
ticities after the anisotropic correction, e*^,^ = (0.791 ± 
3.130, 1.127± 2.331) x 10-^ significantly improves from 
those of raw stellar ellipticities before the correction, 
e,\^ = (-1.179 ±0.014, -0.520 ±0.014) x lO^^. Here, 
the number of stellar objects is iV = 8929. 

Next, we correct the isotropic smearing effect of galaxy 
images to estimate the reduced distortion signal, ga, 



Name'' 



Re 



[mini 



Seeing'^ 
[arcsec] 



' h 



9a — (Pg )ape.p, 



(2) 



comalO 


24.5 


13.75 


0.93 


0.48 


comall 


24.5 


13.75 


0.65 


0.34 


comal2 


24.5 


13.75 


0.63 


0.32 


comalS 


24.5 


13.75 


0.63 


0.31 


comal4 


24.5 


13.75 


0.65 


0.33 


coma20 


24.5 


14.75 


0.75 


0.39 


coma21 


24.5 


13.75 


0.59 


0.29 


coma22 


24.5 


13.75 


0.75 


0.39 


coma23 


24.5 


13.75 


0.63 


0.33 


coma24 


24.5 


13.75 


0.57 


0.28 


comaSO'^ 


24.5 


13.92 


0.69 


0.36 


comaSQf 


24.5 


13.92 


0.67 


0.35 


coma31 


24.5 


13.75 


0.61 


0.30 


coma32 


24.5 


13.75 


0.71 


0.37 


comaSS 


24.5 


13.75 


0.83 


0.54 


coma34 


24.5 


13.75 


0.70 


0.36 


coma41 


24.5 


14.58 


0.72 


0.38 


coma42 


24.5 


14.58 


0.72 


0.38 


coma43 


24.5 


13.75 


0.76 


0.40 


coreS 


42.0 


- 


0.81 


0.41 


sub^ 


16.0 


- 


0.83 


0.41 



Note. — 

'^ data name 

^ exposure times in the R^ and V , respectively. 
'^ the seeing FWHM in unit of arcseconds, for Re band 
the median stellar half-light radius in unit of arcseconds, for iJc 
band. 
'' Re band data taken in March 1st, 2011. 



f 



s Data retrieved from SMOKA. 



Re band data taken in March 30th, 2011. 



where P^o is the pre-seeing shear polarizability tensor. 
The measurement of P^o is very noisy for individual 

faint galaxies because of its nonlinearity (Bartelmann 
& Schneider 2001), which may give a systematic bias 
in weak-lensing distortion measurements. We therefore 
calibrate P^o with the following procedures, in a similar 
way as Umetsu et al. (2010) and Oguri et al. (2012). The 
polarizability tensor is first computed as a scalar polariz- 
ability, {Pg)ai3 = |tr[Pg](5„^. We then compute a median 
for {Pg)ap in Tg, with an adaptive grid to assemble as uni- 
formly as possible. Here, rg is the Gaussian smoothing 
radius used in the KSB method. A sample of galaxies 
satisfies the following conditions to suppress the noise; a 
detection significance level of i^ > 30, a size condition of 
fh > ^h+o'(^h) ^^"^ '^s -* ^5+'''(''g) ^^'^ a positive raw Pg. 
Here, f^ (c(?',*)) and f* (CT(r*)) are the median (rms dis- 
persion) of half-light radii and Gaussian smoothing radii 
for the stars selected above, respectively. We interpo- 
late the polarizability tensor for individual galaxies as a 
function of Vg . The similar interpolation for the absolute 
value of the ellipticity, |e|, is also applied. We use galax- 
ies for shear catalog with :^ > 10 and the same size cut 
as in the calibration. An rms error of the shear estimate, 
ag, is computed from 50 neighbors in the magnitude- rg 
plane. We also assign the wight function for individual 

^ http://www.ifa.hawaii/kaiser/IMCAT 
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objects. 



(3) 



where a is the softening constant variance. We choose 
a = 0.4, which is a typical value of the mean rms ag 
over the background sample (Hamana et al. 2003; Okabe 
et al. 2010b; Umetsu et al. 2010; Oguri et al. 2012). 

To check shear calibration, we use a number of re- 
alistic images of which field-of-view is comparable to 
Subaru/Suprime-cam, kindly provided by Oguri, M.. 
The mock images are generated with different seeing sizes 
(0.51.1 arcsec) and the Moffat profile with power slopes 
3 < /? < 12, using the software GLAFIC (Oguri 2010), 
as described in Oguri et al. (2012). Wc found that a 
multiplicative calibration bias, m and an additive resid- 
ual shear offset, c in Heymans et al. (2006) and (Massey 
et al. 2007) are \m\ < 0.03 and |c| < 2 x lO"'', respec- 
tively, for our typical seeing ~ 0'.'7. 

We combine shear catalog constructed from individ- 
ual images. For the overlapping regions, since the same 
galaxies are detected in different imaging data, we esti- 
mate weighted averages of their position and shear with 
Wg. We compared reduced shear for 1.8 x 10^ overlap- 
ping galaxies and confirmed that the deviation, Aga = 
(-2.16 X 10^6 _t g g X 10-4^ 1.93 X 10-6 ± 6.6 x lO"''), is 
negligible. So far, the number of background galaxies is 
-6.7x 10^ 

2.3. Photometry and Background Selection 

We make a secure selection of background galaxies in 
color-magnitude plane, because a contamination of un- 
lensed member or foreground galaxies in shear catalog 
dilutes weak-lcnsing signals and then leads to the un- 
derestimation of gravitational lensing mass, mainly for 
the central regions (Broadhurst et al. 2005; Okabe et al. 
2010b). 

Photometric catalogs were constructed from the mo- 
saic images using SExtractor (Berlin & Arnouts 1996). 
Our SExtractor parameters is optimized for faint galax- 
ies for shape measurements. We compute the total 
magnitude for each object in the AB-magnitude system 
using the MAG_AUTO parameter and color using the 
MAG_APER parameters, respectively. In color mea- 
surement, we degrade the seeing to the worse image. 
The aperture diameter for MAG_APER is 1.5 times of 
the seeing FWHM. The overlapping galaxies serve as a 
monitor of the offset in magnitude. We introduce an 
additional parameter in each data field to describe an 
offset in two magnitudes, and calibrate with simultane- 
ously fitting bright objects of which magnitudes are less 
than 22 mag. The measurement scatter for faint galaxies 
{Re > 24 mag) is typically less than 0.1 mag. The magni- 
tudes and colors for objects are estimated with weighting 
measurement errors. We then match shear and SExtrac- 
tor catalogs. 

The red-sequence of member galaxies is fitted with a 
linear function, using luminous galaxies (i?c < 18 mag). 
We then define background galaxies whose colors are red- 
der than the red-sequence in the magnitude range of 
20 < i?c < 26 mag (Figure 2). The number of back- 
ground galaxies is reduced to ~ 6 x 10^ after the color 
cut, but it is still by 30 — 60 times higher than those of 




Fig. 2. — Color-magnitude diagram. The color shows the number 
of background galaxies in each pixel (0.1 x 0.1). The white solid 
line denote the red-sequence of member galaxies, fitted by a linear 
function for the bright galaxies. 

clusters at z ^ 0.2, obtained by previous studies using 
two path-band filters (Okabe et al. 2010b). The num- 
ber density, Ug ~ 41.3 arcmin" , is also two or three 
times higher than those for clusters at z ^ 0.2 and our 
previous analysis of Coma cluster (Okabe et al. 2010a, 
Ug ~ 23 arcmin" ). Thus, we can use a correspondingly 
large number of background source galaxies in the case 
of nearby cluster weak-lensing analysis. This is because 
of the following two reasons. First, since colors of red- 
sequence galaxies in clusters becomes bluer with decreas- 
ing their redshifts, the number of galaxies, of which colors 
are redder than those of member galaxies, increases. Sec- 
ond, a large covering area of the nearby cluster increases 
background galaxies. Even if member galaxies are con- 
tained in our background catalog, the dilution effect in 
lensing signals is negligible because the ratio of thou- 
sand members to sub-million backgrounds is tiny. Thus, 
weak-lensing analysis of the nearby cluster has a great 
advantage to compensate low lensing efficiency. 

2.4. Mean Lensing Depth 

Since redshifts of individual galaxies in our shear cat- 
alog are not available, we instead estimate mean source 
rcdshift based on a statistical approach. The lensing sig- 
nal depends on the source redshifts through the distance 
ratio. As a reference, we make use the COSMOS pho- 
tometric redshift catalog (Ilbert et al. 2009) estimated 
by combining 30 broad, intermediate and narrow bands. 
Since the Re band is not available in COSMOS catalog, 
we convert magnitudes from R^ to i' bands. We take into 
account for the filter sensitivities of Suprime-Cam, and 
compute magnitudes in i' band. The probability func- 
tion of redshift, dPwh/dz, for our background galaxies 
selected by the color-magnitude plane (Sec. 2.3) is com- 
puted from matching with COSMOS photometric red- 
shift, with a statistical weight of Wg. The mean distance 
ratio is given by 



(Dis/D,) 



dzdPwh/dzDis/Ds, 



(4) 



where Dg and Dig are the angular diameter distance 
to the sources and between the cluster (lens) and the 
sources, respectively. We obtain {Dis/Ds) = 0.9554. The 
mean source redshift, (z^) = 0.61, is slightly lower than 
those ((z,,) - 0.7 - 0.8) for clusters at z - 0.2 (Ok- 
abe ct al. 2013), because we include many background 
galaxies at lower redshifts. Wc mention that the mass 
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estimate for the nearby system do not strongly depend 
on the redshift distribution of background sources. 

2.5. Luminous Member Galaxies 

We define luminous member galaxies identified in spec- 
troscopic observations in order to compare with mass 
properties. We retrieve luminous galaxies, of which 
magnitude are brighter than i' < 18 mag, from SDSS 
DR8 (Eisenstein et al. 2011), in the region of 4 x 4 
square degree centering NGC4874. Furthermore, mem- 
ber galaxies are selected within the redshift range of 
\z — Zc\ < cr„(l -|- Zc)/c and a = 3000 kms~^, where c 
is the light velocity. To complete the catalog of member 
galaxies in the bright end, we also checked redshifts of 
galaxies in NED^. If they are satisfied with the above 
conditions, we added into the catalog. Luminosities of 
individual galaxies are estimated from apparent magni- 
tudes using the k-correction for early-type galaxies, as- 
suming a single redshift of Zc- 

2.6. Model of LSS lensing 

Weak-lensing mass measurement for the cluster at low 
redshift suffers from lensing signals from background 
galaxies between the cluster and the source redshifts, so- 
called uncorrelated LSS lensing. The three-dimensional, 
inhomogeneous mass distribution causes locally strong 
shear pattern, which potentially gives biases in detec- 
tion and mass measurements of localized objects, such 
as subhalos. We therefore quantify the uncorrelated 
LSS lensing effect on each galaxy in shear catalog, us- 
ing the luminosity and photometric redshift retrieved 
from SDSS DR8 (Eisenstein et al. 2011), following Ok- 
abe et al. (2010a). Galaxies are selected by magnitudes 
i' < 24 mag and photometric redshifts between the clus- 
ter, Zph — Zc > Sz ~ (7^(1 + Zc)/c ~ 0.01, and the 
source redshift. where Zph is photometric redshift of each 
galaxy. We here excluded galaxies spectroscopically iden- 
tified as member galaxies (Sec. 2.5). Masses of individual 
galaxies are estimated using galaxy-galaxy lensing result 
using SDSS data Guzik & Seljak (2002). Using the mass- 
to-light ratio in each band {u'g'r'i' z') derived by stacked 
lensing analysis of galaxies, the luminosity is converted 
into the mass. The masses estimated with different bands 
are used to cross-check and calibrate systematic errors 
in the mass-luminosity scaling relation we adopted here. 
Since the uncorrelated LSS lensing is obtained by in- 
tegrating the effect of light deflections due to galaxies 
at different redshifts along the line-of-sight, we use the 
best-fit scaling relations in mass estimates for individual 
objects in order to quantify an average LSS lensing ef- 
fect. The interior mass structure of each halo is assumed 
to be a universal mass profile found in numerical simula- 
tions, which named 'NFW after the initials of proposers 
(Navarro et al. 1996, 1997). The NFW mass model is de- 
scribed by two parameters of the mass and concentration, 
see details in Appendix C. It is well known that there 
exists a correlation between mass and halo concentration 
(Bullock et al. 2001) predicted by hierarchical structure 
formation scenario. We use the mass-concentration rela- 
tion given by the recent numerical simulation based on 
WMAP5 cosmology parameters (Duffy et al. 2008) to 

^ http://ncd.ipac.caltcch.edu/ 



describe the internal structure. The tangential distor- 
tion signals of individual galaxies are computed on all 
source galaxies. We found that the shear estimated in 
the r' band is consistent with that in the z' band, but 
the estimates in the u' g'i' bands are systematically differ- 
ent, as found in Okabe et al. (2010a). The LSS lensing 

model is therefore defined with gi ' = (g^' -I- g^')/2. 
We note that, since we use the galaxy-galaxy lensing re- 
sult, if there exist groups or clusters behind the cluster, 
our model would fail to incorporate those effects. Since 
this would give a bias on mass measurement of the main 
cluster, we consider it in Sections 3.4.1 and 5. 

3. WEAK-LENSING ANALYSIS FOR SUBHALOS 
3.1. Projected Distributions of Mass and Baryons 

We first make maps of lensing convergence field (k(0)), 
luminosity {l{0)) and number density {n{6)) of mem- 
ber galaxies and the model of the LSS lensing signal 
('^lss(^))- In order to identify subhalos in a model inde- 
pendent way, we reconstruct the projected mass distribu- 
tion following Kaiser & Squires (1993) with a Gaussian 
smoothing kernel. The detail of map making is explained 
in Appendix A. We adopt various smoothing scales in 
the range of rgm = 1, . . . , 5 arcmin, stepped by 0.1 for 
1 — 2 arcmin and 0.2 for 2 — 5 arcmin, to optimize for de- 
tection of subhalos which have various mass properties. 
The definitions of subhalos are described in the next sub- 
section. We here present the correlation between mass 
and luminous matters on the projected distribution. 

Figure 3 shows the significance map, v^ defined by 
kJUk.-, with a smoothing scale of FWHM = 4f (rgm = 
2f4), where the reconstruction error, CTk, is calculated 
over local background galaxies (see Appendix A) and a 
typical value in this smoothing scale is Of^ ~ 7.7 x 10^'^. 
Here we did not take into account LSS lensing model. 
The mass distribution in the central region in which two 
cD galaxies NGC 4874 (a,(5)=(194.898°, 27.959°) and 
NGC 4889 (195.034°, 27.977°) exist is elongated in the 
east and west directions. The clumpy structures are 
found everywhere, but anisotropically distributed. In 
particular, the projected distribution of clumpy struc- 
tures is concentrated 30 — 60f southwest of NGC 4874. 
Some clumpy structures are associated with background 
groups in the literature (alphabets in Figure 3 and Table 

2). 

The left and middle panels of Figure 4 show maps of 
luminosity and number density of member galaxies, over- 
laid with the contours of the projected mass distribution. 
The mass and galaxy distribution are clearly correlated 
with each other. The right panel of Figure 4 is the conver- 
gence map of the LSS lensing model. The signal-to-noise 
ratios of LSS lensing map are at most 0.4cr and thus LSS 
lensing effect accounts for small fraction of the observed 
signal. However, the LSS model fails to describe lensing 
signals around some groups behind the cluster (Table 2). 
The reason is likely that our estimation of LSS lensing 
effect is based on galaxy-galaxy lensing and didn't take 
into account of group-scale or cluster-scale structures. 
This issue will be discussed again in the subsection 3.4.1. 

To quantify the correlations appeared in the maps, 
we compute the pixel-to-pixel coefficients between the 
mass maps {k{6)) and the luminosity {l{9)) and density 
{n{9)) maps for member galaxies. The resultant coef- 
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Fig. 3. — The projected mass distribution with a smoothing scale of FWHM= 4.' in a unit of significance level of i^ = k/o-^. We here 
use the shear without taking into account the LSS lensing effect. The contours of significance start at la and steps with a la value. The 
alphabets and numbers denote the name of known background systems (Table 2) and the name for subhalos (Table 3), respectively. 
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Fig. 4. — Luminosity map for member galaxies spcctroscopically identified in the SDSS DR8 and NED catalog {«' < 18;Sec 2.5). Overlaid 
are the contours of the lensing K-field in units of Icr reconstruction error (Figure 3), without taking into account LSS lensing model. Middle 
: density map for member galaxies. Right : mass map of LSS lensing model estimated from galaxy-galaxy lensing using photometric 
redshifts and luminosities for individual galaxies (Sec 2.6). 



TABLE 2 
Known background systems appeared in mass maps 



ID'' 



Name 



^phot 



Reference 



A MaxBCG J195. 08820-1-26. 78870 0.162 

B GMBCG J195.47315-I-26.95810 0.219 

C MaxBCG J195.47907-I-27.16429 0.208 

D GMBCG J195.34791-I-29.07201 0.189 

E MaxBCG J195. 34617-1-29. 18616 0.170 

F NSC J125939-I-290715 0.189 

G GMBCG J193.96542-I-28. 51557 0.257 

H MaxBCG J193.92901-I-28. 76123 0.259 

I SDSSCGB 06685 0.183 

J WHL J125535. 3-1-273104 0.418 



Koester ct al. (2007) 
Hao et al. (2010) 

Koester et al. (2007) 
Hao et al. (2010) 

Koester et al. (2007) 
Gal et al. (2003) 
Hao et al. (2010) 

Koester et al. (2007) 

McConnachie et al. (2009) 

Wen et al. (2009) 



Note. 



:The identification of background systems in Figure 3. .Photometric redshifts. 
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ROSAT X-isy Map 




Fig. 5. — X-ray surface brightness distribution in the 0.1-2.4 keV 
band from ROSAT X-ray satellite. Overlaid are the contours of 
the mass map with FWHM = 8f3 with taking into account LSS 
lensing model. The contour level starts at la and increases in steps 
of Icr. 



ficients for both (kI) and {ku) change from 0.57 ± 0.08 
{7a) to 0.16 ± 0.02 (Scr) with a decrease of spatial reso- 
lution. Here, the errors are estimated by bootstrap re- 
sampling with 200 realizations of k maps, describing that 
noise peaks are accidentally correlated with smoothed lu- 
minous distributions. In short, the correlation between 
mass and member galaxy distributions is at 7 — 8a lev- 
els. We also computed the coefficients between the mass 
map and the LSS lensing map (klss(^)), and found high 
correlations, 0.51 ±0.08 {6a) for r^m = 5.' and 0.28 ±0.02 
(14cr) for Tsm = 1', respectively. The significance level 
is higher with an increase of resolution, indicating that 
LSS lensing signal caused by small background objects 
makes local shear pattern. 

As shown in Figure 5, X-ray surface brightness distri- 
bution of ROSAT X-iay satellite shows elongated X-ray 
distribution in the central region and an excess X-ray 
flux associated with NGC4839 group in the southwest 
direction (Briel et al. 1992; White et al. 1993; Neumann 
et al. 2001). Overlaid is mass contours with smoothing 
scale of FWHM= 8f33 to compare with diffuse emission 
detected with large PSF of the ROSAT. Ahhough dif- 
fuse X-ray emission from NGC4839 group is associated 
with a clumpy mass structure, all mass structures are not 
necessarily correlated with X-ray features. This point is 
discussed in Section 6.6. 



3.2. Selection and Mass Measurements for Subhalos 

We explore, in a model independent way, subhalo can- 
didates by finding peaks in mass maps reconstructed with 
several smoothing scales. As described in Sec 3.1, maps 
of the observed lensing signals are correlated with those 
of LSS lensing model (Figure 4). In order to securely 
identify cluster subhalos and accurately measure their 
masses, it is crucial to minimize the contamination by 
LSS lensing effect. We therefore calibrate the reduced 

shear in an approximate form of pa " = ga — Qa^^ 
to eliminate LSS lensing effect along the line-of-sight as 
much as possible. We then repeat mass reconstructions 
using the calibrated shear catalog. 

We first select subhalo candidates of which peaks in 
mass maps are above a threshold. The mass maps are 
represented as the convolution of lensing distortion pat- 
tern of cluster mass distribution with smoothing kernels. 
Therefore, Gaussian smoothing scales used for the mass 
reconstruction vary from 1.' to 5.' in order to optimize 
for the detection of subhalos with various mass proper- 
ties. Here, the pixelized k field slightly changes by re- 
construction kernels, like top-hat or wavelet filters. We 
use a significance level, v = K/a^, for a selection of sub- 
halo candidates, where k and a^ are the dimensionless 
surface mass density and the reconstruction error, respec- 
tively. Since the variance and skewness of v histogram 
in the pixels depend on the smoothing scale, we identify 
subhalo candidates above the threshold of three times 
of the standard deviation. The threshold of significance 
level in the highest resolution corresponds to i' > 3.4. 
We first identify subhalo candidates at various smooth- 
ing scales. Then, two peaks appeared between two dif- 
ferent smoothing scales are matched with the condition 
d < FWHM, where d is a distance between two peaks ap- 
peared in different scales, and FWHM is full width and 
half maximum of the larger smoothing scale. The num- 
ber of subhalo candidates becomes 49 at this stage. We 
note that Okabe et al. (2010a) used a mass map with sin- 
gle smoothing scale (FWHM=:2.') and applied the lower 
threshold. Therefore, two of seven subhalo candidates in 
previous paper (ids are 6 and 8) are below more conser- 
vative threshold of this analysis. 

Since we minimized the LSS lensing contribution by 
applying the galaxy-galaxy lensing model, eight known 
background objects (Table 2) are below the thresholds. 
However, the model is not perfect to describe fully LSS 
lensing effect. The other three peaks associated with the 
known background objects (Table 2) are detected with 
the above conditions. One is the background object T 
and two peaks are around the object 'F' (see Figure 3). 
These objects are likely to be groups because the lens- 
ing signals arc stronger than expected from single galaxy 
luminosity. Furthermore, there is a possibility that back- 
ground groups are accidentally superposed with cluster 
subhalos, giving a systematic bias on mass estimates of 
subhalos. This point is discussed in Sec. 3.4.1. 

Next, we measure model-independent projected masses 
(Clowe et al. 2000, see also Appendix B) for shear- 
selected subhalo candidates. This measurement has sev- 
eral important advantages. First, a large number of 
background galaxies are available, because a projected 
mass within a circular aperture radius is computed by 
integrating source galaxies outside the radius. The mea- 
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sured projected mass is a cumulative function of radius. 
Thus, comparing to a tangential distortion profile which 
averages up a tangential component of all the galaxies 
residing in radial bins, this approach suppresses the ran- 
dom noise relevant to the intrinsic ellipticity. Second, 
since the measurement subtracts the background mass 
density surrounding subhalos, we can exclude a con- 
tribution of main cluster mass distribution to subhalo 
masses. Third, mass density of subhalos are expected to 
be close to zero outside tidal radius, measured aperture 
mass corresponds to subhalo mass itself. We note that, 
if mass density profile follows the universal NFW profile 
(Navarro et al. 1996, 1997) without any truncation radh, 
the aperture mass is higher than the spherical one (Okabe 
et al. 2010b). As expected from tidal destruction, the ra- 
dial profile of the projected mass is saturated outside the 
truncation radii, rj. We measure projected masses for all 
the candidates. Since the smoothing kernel for the mass 
reconstructions gives rise to centroid uncertainties of the 
candidates, we determine central position by choosing 
maximize lensing signals within 8.' x 8' box whose center 
is aligned to the map peak position. For accurate mass 
measurements of subhalos having the variety of sizes, it 
is important to explore truncation radii where projected 
mass profile is saturated. We systematically compute 
projected mass profiles by changing background annu- 
lus and then statistically determine the truncation radii. 
Here, inner radius changes from Of? to 14f5 by stepping 
0.'2 and the width is fixed at 3f. The projected mass M2D 
is computed from saturated values, taking into account 
error covariance matrix. The measurement method is 
detailed in Appendix B. We repeated the same analy- 
sis but for different background widths and found that 
the result does not significantly change. Mass measure- 
ments used a considerably large number of source galax- 
ies (4 X 10'^ — 2 X 10^). The number is comparable or less 
than those for main clusters at z ~ 0.2 (e.g. Okabe et al. 
2010b) of which the background number densities are 
rig ~ 5 — 20 [arcmin"^] . Less massive subhalos which are 
detected inside more massive ones should be excluded in 
order to avoid double-counting subhalos. We count the 
ith subhalo using two conditions of the radius rt^i > rtj 
and the subhalo mass M2B,i > -^20 j {i 7^ j)- The num- 
ber of the candidates is reduced from 49 to 39 by this pro- 
cedure. As mentioned above, the LSS model fails to fully 
explain lensing signals of background systems, especially 
of groups scales. Furthermore, since there is a possibility 
to detect mass structures behind the cluster, we conser- 
vatively select the candidates hosting spectroscopically 
identified member galaxies within their truncation radii 
as the cluster subhalos. In the end, we identify 32 peaks 
as dark matter subhalos. Three candidates are associated 
with the background systems (Table 2). Four candidates 
have no optical counter: the two ones are and the others 
arc located around ^ 70 f of the south-east direction and 
the south-east direction, respectively. 

Wc refer to these 32 subhalos as integer numbers, in the 
order of right ascension. The resulting subhalo masses, 
Af2D, range from - 2 x lO^^ /j-iMq to ~ 5 x lO^^ /j-i Mq 
(Table 3). As shown in Figure 6, the radial profiles 
of the projected mass clearly show saturation at some 
outer radii. The subhalos are widely distributed from 
the northeast to the southwest on the sky (Figure 3). In- 



terestingly, the direction connecting between Coma clus- 
ter and Abell 1367 which are parts of Coma superclus- 
ter (Gregory & Thompson 1978) agrees roughly with 
the subhalo distributions. Several massive subhalos are 
associated with well-known, spectroscopically identified 
groups in the cluster (e.g. Mellier et al. 1988; Adami et al. 
2005). Galaxies or groups associated with subhalos are 
summarized with references in Table 3. The cD galaxies, 
NGC4874 and NGC4889, are associated with subhalos 
'21' and '24', respectively. The mean mass ratio between 
this paper and previous paper for overlapping subhalos 
is {Mncw/Moid} = 1-02 ± 0.54. We also measured the 
projected masses for two subhalos of which peaks are be- 
low the threhold in this analysis. The mean mass ratio is 
(Mncw/Moid) = 0.74 ±0.66. Since the number density of 
background galaxies in previous analysis is about half of 
this analysis, we cannot rule out a possibility that these 
peaks are accidentally above the threhold. 

3.3. Stacked Lensing Analysis for Subhalos 

We next conduct stacked lensing analysis for subhalo 
candidates, which is complementary to the projected 
mass measurement. Power of stacked lensing technique 
is to reduce the random noise due to intrinsic cUipticitics 
by increasing a number of source galaxies. Thus, tan- 
gential profiles even for small and less massive subhalos 
can be computed and their average parameters can be 
determined with lower measurement errors. 

We first divide subhalos into three subsamples based 
on the model-independent projected masses of M2D < 
4.6 X 10^2 X 10^2 /i-iMo, 4.6 x lO^^ h-^MQ< Afzo < 
10^3 /i-IMq and lO^^ h-^MQ< Afao- Here, mass 
thresholds are chosen by a subhalo mass function which 
will be described in Sec 6.1. The number of subhalos are 
21, 8 and 3 from less massive to more massive subhalos, 
respectively. As shown in Figure 7, clear breaks are dis- 
played in the stacked lensing profiles. Outside the breaks, 
the profiles are proportional to 6*"^, which indicates that 
the mass density becomes zero outside. We emphasize 
that such a feature was not found in massive clusters 
(e.g. Okabe et al. 2010b) but was found in stacked lensing 
profile for subhalos in our previous paper (Okabe et al. 
2010a). We did not apply any rescaling procedures on the 
radial bins corresponding to lensing signals, because this 
mass weight scheme gives a bias on mass estimates, as 
described by Okabe et al. (2013). Here, an off-centering 
effect (Yang et al. 2006) in lensing signals from the main 
cluster mass is negligible because the mean projected dis- 
tance from the cluster center are much larger than the 
maximum radius. We fit the stacked lensing profiles with 
NFW, TNFW and TNFWProb models (Appendix C). 
Here, TNFW model is truncated NFW (Takada & Jain 
2003; Okabe et al. 2010a). The TNFW model is an ex- 
treme case of truncation models, the mass density outside 
the truncation radius is zero as described in Appendix 
C. The TNFWProb model is the TNFW model taking 
into account a probability function for truncation radius 
which is assumed to be Gaussian with the mean, (rj), 
and the standard error art ■ Given the function, we mea- 
sure a mean subhalo mass (Mgub)- In the model fitting, 
we propagate the systematic errors by possible back- 
ground structures around subhalos '1' and '32', which 
is described in 3.4.1. As expected from clear breaks, the 
mean tangential profiles are well fitted by the TNFW 
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Fig. 6. — The radial profiles of projected mass for subhalos, showing that the masses are saturated at outer radii. The integer numbers 
at the left-top corner denote subhalo names (see also Figure 3). 
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TABLE 3 
Subhalos' properties 



ID'' 



(RA.Dec)'' 
dcg 






iS/N)M, 



Representative galaxies'' 



It 


(193.885 


27.505) 


15.42 ± 2.79 


8.23 


5.98 


NGC4807 


2 


(194.011 


27.685) 


8.79 ±4.69 


6.39 


3.55 


NGC4816 Group" 


3 


(194.197 


27.763) 


3.71 ±1.08 


7.00 


4.61 


SDSS J125645.42+274638.0 


4 


(194.232 


27.053) 


2.89 ±1.08 


5.63 


3.51 


SDSS J125647.00+270324.9 


5 


(194.298 


27.360) 


5.00 ±2.34 


5.35 


3.86 


2MASX J12571076+2724177 


6 


(194.355 


27.660) 


2.52 ± 1.27 


5.93 


4.45 


G12 Group* 


7 


(194.361 


28.187) 


5.99 ±2.84 


5.42 


3.80 


UGC08071, 2MASX J12572841+2810348 


8 


(194.372 


27.380) 


1.87 ±0.73 


4.94 


3.54 


2MASX J12573148+2723048 


9 


(194.381 


27.493) 


12.11 ±2.52 


9.71 


6.45 


NGC4839 Groupt, G4 Group'', NGC4842, X-ray subhalo''^ 


10 


(194.477 


28.507) 


3.24 ±0.75 


8.77 


3.42 


2MASX J12575392+2829594 


11 


(194.572 


27.745) 


4.13 ±0.85 


6.82 


4.03 


2MASX J12581922+274543 


12 


(194.579 


27.846) 


2.02 ± 0.78 


4.56 


3.87 


SDSS J125818. 20+275054.5 


13 


(194.597 


27.101) 


2.70 ±0.77 


6.37 


3.61 


2MASX J12581552+2705137 


14 


(194.640 


27.623) 


4.51 ±1.27 


7.71 


3.53 


NGC4853 


15 


(194.656 


27.905) 


2.96 ±1.44 


5.05 


6.90 


NGC4839 Group" 


16 


(194.659 


26.738) 


5.03 ±1.06 


7.46 


4.19 


SDSS J125839. 93+264534.2 


17 


(194.718 


27.825) 


3.13 ±0.74 


5.60 


4.94 


G9 Group^ SA 1656-030'' 


18 


(194.732 


27.759) 


6.48 ± 2.03 


5.96 


4.47 


G8 Group'' 


19 


(194.790 


28.288) 


4.66 ±1.26 


7.00 


4.74 


SDSS J125914.99+281503.6 


20 


(194.879 


28.062) 


2.90 ±1.58 


3.90 


4.16 


2MASX J12593141+2802478 


21 


(194.882 


27.936) 


4.29 ±1.06 


6.48 


7.23 


NGC4874(cD),part of Gl Group'', X-ray subhalo 2«« 


22 


(194.895 


28.511) 


4.50 ±1.90 


6.75 


3.54 


2MASX J12594129+2830257 


23 


(194.971 


27.837) 


3.75 ± 1.04 


6.48 


4.26 


J194.9353+27.83393t, SA 1656-054^X-ray subhalo 3"" 


24 


(195.052 


28.005) 


5.20 ±2.40 


4.68 


4.71 


NGC4889(cD),part of Gl Group^ X-ray subhalo lHH 


25 


(195.086 


28.542) 


3.86 ± 0.95 


6.09 


3.93 


2MASX J13002268+2834285 


26 


(195.111 


28.654) 


2.75 ± 0.79 


5.54 


4.43 


SDSS J130037.14+283950.9 


27 


(195.115 


28.080) 


4.28 ± 1.74 


4.53 


6.24 


SDSS J130030.95+280630.2,part of G7 Group'' 


28 


(195.155 


28.331) 


5.70 ±1.68 


6.59 


3.68 


NGC4896 


29 


(195.220 


28.010) 


3.64 ±1.30 


6.93 


4.31 


NGC 4908, NGC 4908 Group 


30 


(195.300 


28.558) 


3.12 ±0.66 


7.72 


4.03 


SDSS J130114.96+283118.3 


31 


(195.325 


27.830) 


2.97 ±1.42 


5.59 


3.41 


G4 Group'', NGC4919 


32tt 


(195.421 


29.054) 


45.95 ±7.57 


11.40 


8.35 


G15 Group'', IC 4088,2MASX J13014399+2859587 



Note. — ^ : the name of subhalos. ^ : the wcak-lensing center in units of deg. '^ : subhalo masses in units of W^^h~^MQ. '^ : the 
signal- to- noise ratio computed with the Mq^ profiles. ° : the maximum signal- to- noise ratio appeared in mass maps (k). : the name of 
representative galaxies or optical groups, t : A possibility of overlapped background structure, WHLJ125535. 3+273104 (Table 2). tt : A 
possibility of overlapped background structure, GMBCG J195.34791+29.07201 (Table 2). " : MeUier et al. (1988) ^ : Adami et al. (2005) 
'' : Conselice & Gallagher (1999) * : Adami et al. (2009) «" lAndrade-Santos et al. (2013) ''^ :Briel et al. (1992) 
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and TNFWProb models (Figure 7). The best-fit subhalo 
masses and truncations are listed in Table 4. The best- 
fit mass and truncation increase with increasing model- 
independent projected masses. If subhalo sample were 
all from fake peaks, we could not find these character- 
istic features. We compute the significance probability, 
Q{vl2, x^ijj/2), that the data gives as poor fit as the ob- 
served value of Xmin t>y chance. The NFW model for the 
lowest- and highest- mass sample is rejected with a signif- 
icance level of Q = 4%. Thus, the NFW model is inade- 
quate to describe the tangential shear profile with breaks. 
The mean ratio between the best-fit subhalo masses and 
the mean projected mass of the subsamples, (M2d), are 
((M2D)/A^sub) = 1-02 ± 0.12 for the TNFW model and 
and ((M2D)/(Msub)) = 1.06 ± 0.15 for the TNFWProb, 
respectively. They are in good agreements. 

Next, we repeat the stacked analysis for four sub- 
samples divided by the projected cluster-centric radii 
(0 - 20f,20 - 40f,40 - 60f, and 60 - 80.'). Since tidal de- 
struction predicts that the truncation radius are statisti- 
cally correlated with three dimensional radius, a stacked 
procedure averages out line-of-sight positions for subha- 
los. It provides us with an information of the dependence 
of mean subhalos' sizes on cluster-centric radius. Fig. 8 
displays the mean tangential profiles with clear breaks. 
The TNFW and TNFWProb models give the better to 
fit the stacked profiles (Table 4). Although the NFW 
model fit is acceptable (Q > 10%), the NFW and TN- 
FWProb are preferred by comparing the goodness-of-fit 
of each model. The mean truncation radius increases as 
the projected radius from the cluster center increases. 
The mean mass ratio are {{M2T)) / Msub) = 0.97 ± 0.11, 
and ((M2D)/(Msub)) = 1-00 ± 0.16 for the TNFW and 
TNFWProb models, respectively. 

3.4. Systematic Errors 

We here assess various systematic errors on our subhalo 
analysis. 

3.4.1. Projection effect 

We here examine the projection effect on lensing mass 
measurement of cluster subhalos, caused by background 
groups accidentally superposed along the line-of-sight. 
Although we quantify LSS modeling based on scaling re- 
lations between mass and luminosities, it fails to fully 
describe massive background structures, like groups or 
clusters. As shown in Figure 3, possible background 
structures, J and D (Table 2), are located within two 
subhalo regions of '1' and '32' (Table 3), respectively. 
Although it is in principle very difficult to discriminate 
between them from the observed lensing signal, a differ- 
ence between expected mass density profiles enables us to 
assess a contribution from background structures in ob- 
served lensing signal. Since interior subhalos are tidally 
destructed by their parent halos, it is expected that its 
mass density profile outside the tidal radius sharply de- 
clines. On the other hand, virialized background groups 
or clusters do not show such a feature as long as there is 
no neighboring massive halo. Indeed, tangential distor- 
tion profiles for individual groups or clusters and stacked 
profiles shows a clear curvature as a characteristic sig- 
nature of the NFW prediction and no evidence of trun- 
cation feature (e.g. Johnston ct al. 2007; Okabe et al. 



2010b, 2013; Umetsu et al. 2011; Oguri et al. 2012; Taylor 
et al. 2012). Thus, a model fitting to tangential distor- 
tion profiles helps us to discriminate between subhalos 
and background objects. As mentioned above, the pro- 
file for less massive subhalos is very noisy and slightly 
changed by a choice of radial bins, because the num- 
ber of source galaxies is small in proportion to the area 
surrounding less massive subhalos. We therefore concen- 
trate on computing distortion profiles for three massive 
subhalos ('1', '9' and '32'), of which masses are higher 
than 10^3 /i-iAfQ. 

Fig. 9 displays breaks in the tangential shear profiles. 
The slope follows ex 9~^ outside the break, as shown in 
stacked lensing analysis (Sec 3.3). The off-centering ef- 
fect (Yang et al. 2006) of main cluster mass on lensing 
signal is negligible because of large separation from the 
cluster center. We first fit TNFW model as a model of 
subhalos to the tangential distortion profiles. We also 
tried to fit a truncated singular isothermal sphere (TSIS; 
Okabe et al. 2010a) model to the profile but found a poor 
fit for massive subhalos. The best-fit masses and trun- 
cation radii are shown in Table 5. We found the best- fit 
values do not change significantly by a choice of radial 
bins. The solid lines for the best-fit values well describe 
the profiles with the breaks (Figure 9). We next fit the 
NFW model to the data and then obtain larger mini- 
mum reduced Xmin values than those for TNFW model. 
The significance probability, Q, for the NFW model are 
0.12, 0.74 and 0.08 for the subhalos '1' , '9' and '32', 
respectively. When we adopt the threshold of 10%, the 
NFW model for the subhalo '32' is unacceptable and for 
the subhalos '1' and '9' are acceptable though Q for the 
subhalo '1' is close to the threshold. If lensing effect 
by background objects accounts mainly for observed sig- 
nals around the subhalos '1' and '32', the profiles should 
be well described by the NFW model. We repeat the 
tangential fits using photometric redshifts 2;phot = 0.418 
and 0.189 as redshifts of background objects around the 
subhalos '1' and '32', and obtain the virial masses for 



r(i) 



the NFW model, M^{> = 12.111^^9^ x IO^/i-IMq and 

M^.p = S.dStlil X 10"/i-iMo, respectively. In this 
mass scale, no clear truncation radius was found in the 
tangential profiles (e.g. Okabe et al. 2010b; Oguri et al. 
2012). Thus it implies that all lensing signals cannot 
be explained solely by background objects. We then si- 
multaneously fit with the TNFW and NFW models for 
subhalos and backgrounds, where we assume the mass 
and concentration relation for backgrounds (Duffy et al. 
2008) and fix the truncation radius derived by fitting 
with the TNFW model. Although measurement errors 
of subhalo masses become larger, we find that the best-fit 
subhalo masses are decreased by 30-40%. We therefore 
add these differences to the second error in Table 5, as 
the systematic error. 

3.4.2. Probability of Spurious Peaks 

To measure a reliable subhalo mass function, it is of 
critical importance to statistically rule out a possibility 
that these subhalo candidates are spurious peaks. The 
peak finding method always suffers from the presence of 
spurious peaks. It is therefore of vital importance to 
quantify the number and properties of spurious peaks 
in order to confirm the purity of real subhalos. For 
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Fig. 7. — The mean distortion profiles obtained by azimuthally averaging the measured galaxy ellipticities for 32 subhalos. The subsamples 
are selected with the model-independent projected masses. Left : 21 subhalos with M20 < 4.6 X 10^^ h~^MQ. Middle : 8 subhalos with 
4.6 X IOI2 /x-lM0< M20 < 4.6 X IO12 h-'^MQ. Right : 3 subhalos with 10 X 10^^ h-'^MQ< M2D. The profile slopes drastically change 
at the truncation radii. The red solid, green dashed and blue dotted curves are the best-fit TNFW, TNFWProb, and NFW models, 
respectively. The NFW and TNFWProb well describe a sharp truncation, while the NFW model for less massive (left) and massive (right) 
subsamples is strongly disfavored. The best-fit truncation radius increases with increasing the mass of the subsample. 

TABLE 4 
Best-fit mass parameters for TNFW and TNFWProb models for lensing-selected subhalos 



Sub-sample^ 




A^Bub 






(A4nb)^ 


in)" 

h-^kpc 


h-^kpc 10 


"'h-^L^ 


,0 


n-ake*^ 


-4.6 X 10^^ h' 
(4.6 - 10] X 10 
10^^ H^^Mq - 


'^- h- 

J 


t 21 

1j\/q 8 
3 


2.9lt^;^2^ 

5.93+|.« 

26.72_4;j^g_5 gg 


27.48+flt 
70 70+25.42 

I6l.l5tg;i 


n c;q+U.4y 

d.Od_Q 44 

5.95ti;«« 

23.13_g;33_5 Qg 


35.57tl-i| 

70 10-1-33. 06 
(0.iO_22 20 

127.1lt^i;f2 


9.97-:- 
< 1.75 

nn 01 +18.39 
00. 0±_j^gg5 


2.11 
5.24 
7.49 


13.69 
8.67 
5.35 


0.019 

0.011 

2 X 10"^ 


- 20.' 
20.' - 40f 
40.' - 60f 
60.' - 80f 




11 
10 
8 
3 


■1-U0_Q gg 

tr 40+1.04 
^•^■'-1.33-0.73 

"^0 9Q+3.21 


35.101^'^^ 
A9.29tf,lt 

209.69+^3''^g 


■^-UO.g g2 

c r.n+0.7:i 
^•^^-1.09-0.17 

30.27+l?ti.75 


35.66ifg^^ 

40 oq+9.77 

49.56tl?:I^ 


< 19.71 

< 19.43 

U.62+_ni 

< 19.00 


5.35 
3.47 
1.52 
6.16 


8.25 
10.70 
8.39 
7.62 


0.082 

0.062 

0.023 

5 X 10-^ 



Note. — ^ : the name of subsamples for subhalos in stacked lensing analysis. : the number of subhalos. '^ : the best-fit subhalo 
mass and truncation radius for TNFW model. '^ : the best-fit subhalo mass, and the average and standard error of truncation radius 
for TNFWProb model. '^ : the average luminosity for associated galaxies, estimated by weighting tangential distortions, g+. ' : the 
signal-to-noise ratio for tangential distortion profile. ^ : the probability that the TNFW mass and truncation radius for fake subhalos are 
within Icr contours for the best-fit values of observed ones, t ; the NFW model is strongly disfavored. 



TABLE 5 
Best-fit masses and truncation radii for three massive subhalos using the tangential distortion profiles 



ID'' 



1 
9 
32 



lO^'^h-'^M, 

+2.37 







14.26 



2.53-5.55 

11.05^fg4 



/i^^kpc 

+2.V4 



x4in/d-o.f^ S/N- 



Xi^jd.o.i (NFW)* 



77.22 

68.70; 

184.38" 



-3.81 

-5.54 

-9.49 

-14.75 

-16.65 



5.38/3 

0.59/3 

3.74/5 



6.17 
5.99 
8.33 



8.79/4 
1.97/4 
11.31/6 



Note. — ^ : the name of subhalos (Table 3). ^ : the best-fit mass in units of 10*^^ h~^MQ. '^ : the best-fit truncation radii in 
units of /i~^kpc. '^ : The reduced chi-square for the best-fit truncated NFW (TNFW) model and d.o.f is the degree of freedom. "^ : The 
signal-to-noise ratio for the tangential distortion profile. * : The reduced chi-square for the best-fit NFW model. All are higher than those 
of TNFW model. In particular, the significance probability, Q, for the NFW model of subhalo '32' is less than 10%, indicating that the 
profile is not well fitted by the NFW model. 
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Fig. 8. — The mean distortion profile for four subsample of 32 subhalos, selected by the projected cluster-centric radii of — 20.' (top-left), 
20 — 40f (top-right), 20 — 40.' (bottom- left) and 60 — 80.' (bottom-right). The lines have the same meaning as those in Fig. 7. 
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Fig. 9. — Radial profiles of tangential shear component (top panel), g_)_, and the 45 degree rotated component (bottom panel), gx, for 
three massive subhalos ('!', '9' and '32'; from left to right). The tangential signals sharply decline outside truncation radii. The red solid 
and blue dotted lines are the best-fit for TNFW and NFW models, respectively. 



this purpose, we create 200 bootstrap data-sets gener- 
ated by randomly swapping reduced shear at fixed po- 
sitions, repeat map makings and then identify artificial 
fake peaks satisfied with the same condition except for 
spectroscopic information of galaxies. The number of 
spurious peaks is 5.32 ± 2.23 per each realization. Thus, 
wc cannot perfectly rule out a fake contamination in the 
subhalo catalog. It is however difficult to quantify the 
purity of subhalos comparing the numbers of detected 
subhalo candidates and fake peaks, because we excluded 
several subhalo candidates taking into account the cloud- 
in-cloud problem and background groups. On the other 
hand, stacked lensing analysis enables us to measure 
mean parameters even for spurious peaks, though indi- 
vidual measurements of spurious peaks are very difficult 
due to high statistical noise. By comparing statistical 
properties of spurious peaks, it allows us to discuss the 
purity of our subhalo catalog. We generated 500 boot- 
strap replications of stacked tangential profiles using the 
catalog of artificially fake peaks. Here, the number of 
spurious peaks and the radial bins are the same as those 
for subsamples in stacking lensing analysis of subhalos 
(Sec 3.3). The mean mass and truncation radius for 
spurious peaks are estimated by fitting with the TNFW 
model. We found that 38 — 45% of profiles for spuri- 
ous peaks gives a poor fit {Q < 10%) or are ill con- 
strained. It indicates that the profile shapes are different 
from those of observed subhalos. Indeed, the stacked 
lensing analysis for shear-selected subhalos shows that 
the best-fit truncation radius depends on the projected 
mass and the cluster-centric radius (Figures 7 and 8). 
If our sample were all spurious peaks, such a clear de- 
pendence could not be found. To make a more robust 
conclusion, we estimate a probability that these parame- 
ters for spurious peaks accidentally coincides with those 
for observed subhalos within 1 a uncertainty, based on 
Monte Carlo re-distributions of the best-fit values with 



the error covariance matrix. The fake probabilities for 
individual subsamples in stacked lensing analysis , T'fake , 
are from 10-^% to ~ 8% (Table 4). Muhiplying the 
number of subhalos by the fake probability in each sub- 
sample, the expected numbers of spurious subhalos are 
less than unity. Therefore, we conclude that our sample 
of subhalos has high purity. 

3.4.3. Selection Criteria 

A choice of the threshold in the signal-to-noise ratio 
for mass maps gives one of main systematic errors, be- 
cause we cannot completely rule out a possibility that 
peak heights in mass maps are accidentally above or be- 
low the threshold due to reconstruction errors of finite 
sampling of background galaxies. Assuming Poisson fiuc- 
tuations of the reconstruction noise, the signal-to-noise 

ratio changes by S{S/N) = ^.JTVbkg/^bkg = N'^^'^ ~ 
0.25, where iVbkg is an effective number density in the 
smoothed mass map with the highest resolution. We re- 
peated the same analysis with different thresholds. The 
sample numbers become 24 and 49 with higher and lower 
thresholds, respectively. This systematic error is taken 
into consideration to compute a subhalo mass function 
in Sec 6.1. 

3.4.4. Stacking Procedure 

We investigate whether the stacking method gives 
systematic errors, because the mean tangential profile 
stacked over subhalos with various truncation radii would 
make the break feature blunt. We make synthetic weak 
shear catalogs of subhalos with analytic TNFW model 
and the intrinsic ellipticity. Here, the number of back- 
ground sources are the same as our observation. The 
parameters of TNFW model for individual subhalos are 
generated from a Gaussian distribution. The mean and 
standard error of subhalo masses and truncation radii for 
simulated samples are shown in Table 6. We assume that 
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the coefScient between subhalo mass and truncation radii 
is 0.7 and the halo concentration is 1, for the simphcity. 
We compute 500 samples of each stacked profile and fit 
with the TNFW and TNFWProb models. As shown in 
Table 6, the mean tangential profiles recover well the in- 
put values. In TNFWProb model fitting, since the mean 
and standard error of the truncation radii are degener- 
ate, (rt) cannot be well constrained in some cases and 
then the mean measurement error become large. 

4. CLUSTER GALAXY - GALAXY LENSING ANALYSIS 

Galaxy-galaxy lensing analysis for member galaxies se- 
lected solely by their luminosities provides us with com- 
plementary and important information of cluster subha- 
los, because the sample is unbiased with respect to the 
lensing definition of subhalos (Sec 3). We here use mem- 
ber galaxies of which luminosities in the z' band are larger 
than 10^'^ h~^ Li' ^Q. We compute a stacked tangential 
distortion profile as a function of a radius from lumi- 
nous member galaxies. Since the cluster field is crowded, 
neighboring luminous galaxies would be a serious con- 
tamination in the stacked lensing profile, if they are not 
separated enough. We thus need to determine an outer- 
most radius of the profile in order to minimize the lens- 
ing contamination from the neighboring luminous mem- 
bers. We estimate the histogram of projected distances 
between the luminous member galaxies. The outermost 
radius is chosen to be 5' by applying a threshold that 
the mean number of neighboring luminous galaxies is 
less than the unity. We compile 64 luminous member 
galaxies located in the projected cluster-centric radius of 
lOf < r < 80f. Figure 10 shows the mean tangential pro- 
file. A sharply truncated profile is not found in contrast 
to lensing-selected subhalos (Figures 7 and 8). The NFW 
and TNFWProb models are thus applied to describe the 
profile. Two models (Table 7) give an acceptable fit. 
The best-fit tangential profile for the TNFWProb model 
is similar to that for the NFW one, because a standard 
error of the rt distribution makes a sharply truncated 
profile blunt. The mean virial radius for the NFW model, 
273.87t45'98 'i~"'^kpc, is larger than the mean truncation 

radius (r*) = 60.84^1^4 h'^kpc for the TNFWProb 
one. Three luminous galaxies are on average inside the 
mean virial radius of luminous galaxies. In other words, 
the mass distribution of subhalos associated with lumi- 
nous galaxies are overlapped with each other. The NFW 
model is therefore unlikely to represent halos associated 
with luminous member galaxies. On the other hand, the 
TNFWProb model gives a large scatter of the trunca- 
tion radius compared to the mean value, art/{rt) '^^ 37%. 
A broad distribution of the truncation radius smoothcn 
truncation feature in the mean tangential profile for a 
large sample of less massive subhalos, which makes us 
difficult to resolve the subhalo size. To confirm this ex- 
planation, we conducted the stacked lensing analysis us- 
ing mock shear catalog (Table 6) in the same way as Sec. 
3.4.4. We found that truncation feature in stacked profile 
is obscured and the mock simulation recovers the input 
values. Although we also conducted fitting stacked pro- 
files for subsamples divided by luminosities or projected 
distances from the cluster center, only upper limits can 
be derived. We also checked less luminous galaxies of 
which luminosities are less than lO^"/i"^X0, but found a 
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Fig. 10. — The mean distortion profile for 64 luminous member 
:alaxics of which luminosities arc more tha 
have the same meaning as those in Fig. 7. 



galaxies of which luminosities are more than 10^'^ /i ^Lq. The lines 



serious contamination from neighboring luminous or less 
luminous galaxies in the mean tangential profile. 

5. MAIN CLUSTER MASS MEASUREMENT 

A tangential distortion profile, g+, with respect to the 
cluster center is powerful to estimate the cluster mass 
(Okabe et al. 2010b). The tangential distortion profile 
as a function of the projected cluster-centric radius is 
computed by azimuthally averaging the measured galaxy 
ellipticities. Here, we assume that the cluster center is 
the central position of the subhalo '21' associated with 
the cD galaxy (NGC4874). Top panel of Fig. 11 shows a 
complex feature of the lensing profile in the radial range 
of If — lOOf, extending over two orders. Here, the data 
points are calculated by shear catalog without the LSS 
modeling. The lensing signal changes from 0(10^^) to 
0(10^"^) as the radius increases. As expected from low 
lensing efficiency of the nearby cluster, the lensing sig- 
nal is one order lower than those for massive clusters at 
z ^ 0.2 (e.g. Okabe et al. 2010b). However, the signal-to- 
noise ratio reaches S/N ~ 13.3 thanks to a remarkably 
large number of of background galaxies, which is compa- 
rable to or higher than those of clusters at ^ ~ 0.2 (Okabe 
et al. 2010b). This high signal-to- noise ratio certifies a 
weak-lensing analysis for low redshift clusters (z ^ 0.1) 
which have been overlooked for a long time. We also find 
that the 45 degree rotated component, gx (bottom panel 
of Fig. 11), which is non-lensing mode serving as a mon- 
itor of systematics errors, is one order smaller than the 
lensing mode, g+. It is consistent with the null signal. 

The tangential distortion contains the full information 
on the lensing signals, such as from the smoothed mass 
component of the main cluster, the interior substructure 
(Sec 3.2) and LSS lensing signals behind the cluster (Sec 
2.6). To understand the profile, we compute the tangen- 
tial shear profiles for 32 shear-selected subhalos and LSS 
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TABLE 6 

Results of stacked lensing analysis using mock shear catalog 



A^s, 



(A^uD'^ 



^(A^input^b 



sub 



\b / input \ b / input \ 



nput-,b 



21 


3.5 


8 


7.0 


3 


20.0 


11 


3.0 


10 


5.0 


8 


5.0 


3 


30.0 


64 


1.0 



1.0 
2.0 
5.0 
1.0 
1.0 
1.0 
5.0 
0.5 



35 
70 
130 
35 
50 
50 
200 
60 



10 
20 
40 
20 
20 
15 
20 
20 



1.01 ±0.09 (0.11) 
0.92 ±0.21 (0.21) 

1.02 ±0.26 (0.21) 
0.98 ±0.28 (0.17) 
0.92 ±0.17 (0.18) 
0.93 ±0.17 (0.19) 
1.04 ±0.26 (0.24) 
0.91 ±0.27 (0.47) 



1.10 ±0.12 
0.96 ±0.36 
1.08 ± 0.49 
1.18 ±0.38 
0.99 ±0.33 
1.00 ±0.35 
1.10 ±0.28 
1.00 ±0.20 



(0.21) 
(0.29) 
(0.35) 
(0.79) 
(0.37) 
(0.36) 
(0.52) 
(0.77) 



1.03 
1.04 
1.03 
1.02 
1.01 
1.02 
1.04 
1.18 



±0.12 (0.15) 
±0.23 (0.23) 
±0.16 (0.20) 
±0.18 (0.22) 
±0.19 (0.21) 
±0.18 (0.22) 
±0.14 (0.15) 
± 1.18 (1.00) 



1.12 ±0.18 (0.29) 
1.11 ±0.44 (8) 

1.02 ±0.21 (0.6) 
1.11 ±0.40 (5) 

1.09 ±0.35 (13) 
1.11 ±0.41 (5) 

1.03 ±0.11 (11) 

1.04 ±0.75 (6) 



Note. — ^ : the number of simulated subhalos. ^ : the mean and standard error of the mass (10^^ 
(/i~^kpc) for simulated subhalos. '^ : the mean ratio of the outputs to the inputs for TNFW model. The 
deviation based on 500 realizations. The values in the brackets are the mean measurement uncertainties. '^ 
to the inputs for TNFWProb model. 



h~^MQ) and truncation radius 
errors present are the standard 
: the mean ratio of the outputs 



TABLE 7 
Best-fit mass profile parameters for NEW and TNFWProb models, obtained by galaxy-galaxy lensing for luminous 

MEMBER galaxies 



Number^' {M^h)^ (rt)^ 

IQl^fe-^AfQ fe-lkpc 



64 



1.10 



-U.4U 
-0.40 



60.84 



+'IM1 — 
15.34 



fe-lkpc 10l2/i-lMfr 



22.56 



+ia.*/ 

10.64 



2.39 



+ 1.89 
1.01 



22.99 



9.81 






S/N 



2.25 



5.63 



Note. — ^: the number of luminous member galaxies selected by luminosities (L^/ > 10^'^/i^^Lq) and the cluster-centric radii 
(10.' < r < 80.'). ^ : the best-fit mass, the average and the standard error of truncation radius distribution for TNFWProb model. '^ : the 
best-fit virial mass and halo concentration for NEW model. '^ : the average luminosity {L^i). '^ : the signal-to- noise ratio for the tangential 
distortion profile. 



lensing model, described by blue dotted and magenta 
dashed-dotted lines in top and middle panels of Fig. 11. 
The signal-to-noise ratio for the subhalos is S/N ~ 4.4, 
accounting for 33% of the total distortion signal. It indi- 
cates that the profile highly resolves lensing signal from 
the interior substructure by the large apparent size of the 
cluster. The observed signal in a central region (r ^ b'.) 
is dominated by the subhalo '21'. The lensing signals 



in 5.' 



12 f and in r ~ 70f arc depressed by promi- 



nent subhalos. As for LSS lensing, the signal-to-noise 
ratio is S/N ~ 1.3. Here, since the LSS lens modeling 
for the possible background groups, 'F' and T (Table 2 
and Fig. 3), is significantly failed, we estimated the lens- 
ing distortion pattern from the NFW profile determined 
by tangential profile for this object and found that this 
background group depresses the observed lensing signal 
at r ~ 70f. 

We fit single NFW model to the distribution profile 
corrected with the LSS lensing model. The lensing sig- 
nals from the subhalos gradually changes from positive 
to negative in the central region ( ^ 12f). These abso- 
lute values account for a large fraction of the total lensing 
signals, which makes it difficult to discriminate between 
signals from the subhalos and the main cluster from the 
central signals. Indeed, when a single NFW model fits 
to our data, the best-fit halo concentration is systemat- 
ically changed (Acvir ^1) by a choice of the innermost 
radius. To avoid the subhalo bias in cluster mass mea- 
surement, we estimate the radial range for the fitting, 
by applying the fraction of the absolute value of subha- 
los signals to the observed ones is less than 30%. We 



fit the tangential profile between 13f (~ 260 h ^kpc) 
and 64f (~ 1.3 /j-^Mpc) with single NFW model. Here, 
the physical scale of the innermost radius is compara- 
ble to those for massive clusters aX z ^ 0.2 (Okabe 
et al. 2010b). The resultant masses at different over- 
densities are listed in Table 8. The virial mass and con- 
centration are M^i,= d,A2t2A2 ^ lO^'^h^^ h~^MQ and 

Cvir= 3.57j;};^2, respectively. 

We also conducted fitting the profile in the full range 
with the NFW model for the smooth matter component 
and a central point mass contribution of the BCG. We 
obtain the point mass Mpt = 5.67to;32 x 10^^ H'^Mq, 
which is consistent with the projected mass measurement 
of the subhalo '21' (Table 3). The summation of the 
viral mass and the point mass is Mtot = 8.2llj^gg x 
10^** h~^AlQ. The total mass is in a agreement with the 
viral mass estimated derived by the tangential fit for the 
radial range to minimize subhalo contributions. 

We next repeat fitting with the NFW model as the 
smooth component, by fully taking into account lens- 
ing signals from all shear-selected subhalos and LSS 
lens model. The best-fit profile for the smooth compo- 
nent is show in green dashed line of Fig.ll. The to- 
tal signal (red solid line) from three different compo- 
nents of the smooth NFW profile, subhalos and LSS 
lens model describes remarkably well the observed sig- 
nals. The summation of the viral mass and subhalos' 
masses, Mtot = 8.18;2;o2 x lO^** h~'^MQ, are in a good 
agreement with the virial mass M^ir (Table 8). A singu- 
lar isothermal sphere (SIS) model is strongly disfavored 
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as the smoothed mass component, returning a goodness- 
of-fit statistic of Xmin/d-o.f. = 59.3/11. 

The projected mass {M(^J measurement (Cc statis- 
tics) is less sensitive to subhalos, because it estimates 
a cumulative profile. It is thus complimentary to the 
tangential fit. Figure 12 shows the M^^ profile calcu- 
lated with the fixed background annulus of 70 f — 90 f. 
The background region is inside the best-fit virial ra- 
dius (rvir = 96.67tio;32 arcmin) derived by the tangen- 
tial shear fit. Following Okabe & Umetsu (2008), we fit 
the C,c profile with a single NFW model, taking into ac- 
count the error covariance matrix. The best-fit profile is 
shown by red solid line. The best-fit values (Table 8) are 
consistent with those for the tangential shear fit. 

We compare the best-fit mass and concentration with 
the literature (Table 9). Our mass measurements are 
consistent with our previous analysis (Okabe et al. 
2010a). A statistical precision of new mass estimates 
improved by four times thanks to a huge number of 
background galaxies. Gavazzi et al. (2009) conducted 
fitting the tangential distortion profile in the range of 



r [h ^Mpc] 



0f35, 



35.' with the NFW model including and ex- 



cluding priors on mass and concentration relation, us- 
ing CFHT/Megacam data. They didn't take into ac- 
count the LSS lensing effect. The best-fits, regardless 
of priors, are compatible within their large errors with 
ours. However, their lensing signals (Fig 3 in Gavazzi 
et al. 2009) differs at aU from ours (Fig 11). Their 
profile (0.'35 ^ r ^ 35.') is well expressed by single NFW 
model. Our lensing signals in the same radial range are 
dominated by prominent subhalos in r ^ 12.' and by the 
smooth mass component in 12 ^ r ^ 35f, respectively. 
We conduct the tangential fit using the profile computed 
with the same radial bins as Gavazzi et al. (2009) and 
single NFW model and obtain only an upper limit on 
Mvii< 4 X 10^^ h~^MQ because of an inadequate model. 
Kubo et al. (2007) carried out weak-lensing analysis us- 
ing SDSS data. The best-fit M200 (Kubo et al. 2007) 
derived by fitting the tangential profile up to 10 /i~^Mpc 
with a single NFW model is three times higher than ours. 
They didn't consider background LSS lensing effect, ci- 
ther. Since their outermost radius (10 h~^Mpc) is five 
times higher than our best-fit virial radius, their mass 
would be overestimated by mass distribution outside the 
cluster. The dynamical mass estimates of M^ir and M200 
(Rines et al. 2003; Lokas & Mamon 2003) agree with our 
best-fits, albeit their concentration is three times higher 
than our estimates. 

6. DISCUSSION 

6.1. Subhalo Mass function 

Subhalo mass function is computed from Monte Carlo 
redistributions of subhalo masses taking into account 
both the measurement uncertainties and the systematic 
error on the threshold we applied (Sec 3.4.3). Figure 
13 shows the subhalo mass function. The error bars 
are based on both measurement and systematic errors. 
The resulting subhalo mass function covers over two or- 
ders of magnitude in mass. The number of subhalos de- 
creases as their masses increases, while it is significantly 
dropped in the low-mass end (Msub/-^vir ^ 10^^) be- 
cause of the detection limit appeared in the mass map. 
To access the purity, we also compute the mass function 
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Fig. 11. — Top panel: the tangential distortion component, g-^., 
with respect the projected cluster-centric radius, in the range of 
1 — 100 arcmin, estimated by azimuthally averaging the measured 
galaxy ellipticities. The green dashed, blue dotted and magenta 
dashed-dotted lines are the best-fit NFW profile as the smooth 
mass component of the main cluster, lensing signals expected from 
subhalos and LSS lensing model, respectively. The total lensing 
signals (red solid line) of three components well describes the ob- 
served distortion profile. Middle panel: the tangential profiles for 
subhalos (blue dotted) and LSS lensing model (magenta dashed- 
dotted), are the same as the top panel except in a linear scale. 
Bottom panel: the 45 degree rotated component, gx, is consistent 
with a null signal. 
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Fig. 12. — M(^^ profile as a function of the projected cluster- 
centric radius, estimated by azimuthally averaging the galaxy el- 
lipticities with a correction of LSS lens model. The red solid line 
is the best-fit NFW profile. 
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TABLE 8 
Main cluster mass estimates 



Fitting Method'' 






A/20(J 

10" h-^Mp, 



C200 



+4.17 
2.42 

+2.42 
1.82 



+ i.b4 
-1.12 
+0.80 



+2.53 
1.58 
1.51 
1.20 



+ 1.17 
-0.84 
+0.61 



A/500 

10" H-^Mq 



+ 1.U4 
-0.76 
+0.69 



A/loOO Af2500 

10" h-^MQ 10" h~^MQ 



2 47+U'44 
^•^'-0.37 
Q Q7+0.36 
^•^' -0.33 



, it;+U.22 
^•^'-'-0.22 



3+ profile 
C,c profile 



8.42 
8.31 



3.57 
3.24 



6.23 
6.08 



2.55 
2.30 



3.89 
3.67; 



Note. 



profiles for the fitting 



TABLE 9 

Mass and concentration in the literature 



Reference 


A/vir 

10^'' h-^MQ 


Cvir 


Af200 

10^'' h-^MQ 


C200 


WL : (Kubo et al. 2007) 
WL : w/o priors'* (Gavazzi et al. 2009) 
WL : w/ priors'* (Gavazzi et al. 2009) 
WL : (Okabe et al. 2010a) 
Dynamics : (Rines et al. 2003) 
Dynamics ; (Lokas & Mamon 2003) 


4.27t»-l^ 
7.77+113^ 
8.92^26^?^ 

8.45 ±3.15 


6-7+^3 
9.4 


18.8^5-.;^ 

t;7+3.01 
■^■^'-1.47 
fi 7Q+4.27 

pf;.+12.06 

7.85 


04+13. lb 

■3-»"t-1.84 
q c + 1.1 

2.5011:^1 



Note. — '* priors with and without on the mass-concentration relation. 



of spurious peaks, dA'^fake/dluMfake- Here, we calculate 
the mass function for spurious peaks, using the proba- 
bility distribution of best-fit masses, Mfakc, derived by 
stacked Icnsing analysis (Sec 3.4.2). The number of fake 
peaks as a function of the mass is given by 7V(A'/fakc) — 

./Vfake Y.i ^samplo,i^j(Mfakc)/ Y.i -^sample, i, where TVfake is 

the total number of spurious peaks and Pi(Mfako) and 
A'sampic.i arc the probability distribution and the sam- 
ple number for the subsample divided by the projected 
masses, respectively. The green dashed lines show a sin- 
gle peak of the mass function for spurious peaks. The 
functional form is different at all from the observed mass 
function. The peak hight of spurious peaks is one or- 
der lower than the observed mass function at the same 
mass range. Stacked lensing analysis of fake peaks (Sec 
3.4.2) disfavors the contamination of spurious peaks in 
our sample of subhalos. Even if they exist, the contami- 
nation level is negligible for a study of the mass function. 
We fit subhalo mass function with single power law 
model (e.g. Gao et al. 2012), 



dn/dlnA4ub«M-^ 



(5) 



and a Schechter function (e.g. Schechter 1976; Shaw et al. 
2006), 



dn/dlnMsnh oc M^^^ exp(-Msub/A^») 



(6) 



The cutoff mass, M,, in the Schechter function is sensi- 
tive to the abundance in high-mass end. However, since 
the abundance of massive subhalos is small, it is not well 
constrained, M*/Mvir = 0.099to:TO8- We are therefore 
unable to discriminate between the single power law and 
the Schechter function. The best-fit power indices, which 
characterize the shape of the function for intermediate 
and low range, are in good agreements (a = l.ll^Qjg 

and j3 = 0.94J^Q'2g). We also computed subhalo mass 
function including four subhalo candidates without no 
optical counter and obtain the best-fit a — l.\T^^^ 



and /3 = 1.02to;32. For further check, we excluded the 
most and least massive massive subhalos to construct 
mass function and found that best-fit slopes do not sig- 
nificantly change. The best-fit slopes are in remarkable 
agreements with CDM predictions ^ 0.9 — 1.0 from nu- 
merical simulations (e.g. Diemand et al. 2004; De Lucia 
et al. 2004; Gao et al. 2004b; Shaw et al. 2006; Angulo 
et al. 2009; Giocoh et al. 2010; Klypin et al. 2011; Gao 
et al. 2012) and analytical models (e.g. Taylor & Babul 
2004; Oguri & Lee 2004; van den Bosch et al. 2005; Gio- 
coli et al. 2008). Recent high- resolution numerical sim- 
ulation (Gao et al. 2012) found that the slope of mass 
function in the range of 10^^ < i\/sub A^20o < 10^'^ gives 
a = 0.98. 
The mass fraction for observed subhalos is estimated as 



/sub = T,^ Msub,i/A^v 



0.226 



+0.111 
-0.085 



with the tangential 



fit for the main cluster and fsnh = 0.229lQQg4 with (c 
fit, respectively. Shaw et al. (2006) estimated the mean 
mass fraction as a function of the virial mass, (/sub) = 

0.14± 0.02(Afvir /8 X 10" /i-1Mq)°-44±006^ rpj^g j^g^gg 

fraction for the Coma cluster is larger than the mean 
fraction, (/sub) = 0.16, estimated with our best-fit virial 
mass. We also calculated the mass fraction within r2oo, 
/sub.2oo = 0.222 ±0.077. 

6.2. Correlation Between Subhalo Masses and 
Truncation Radii 

It is interesting to investigate the correlation between 
the subhalos' masses and truncation radii, because they 
are free parameters in stacked lensing analysis. We com- 
pile the stacked lensing results divided by mass and 
cluster-centric bins and find a tight correlation of Mguh oc 

rl^^-'>-°^ and A/,ub oc rl'^^-°-^'<> for the TNFW and TN- 
FWProb models, respectively. Considering a functional 
form of NFW model, Mnfw(< x) oc log(l+a;)-a;/(l+a;) 
(cq. C2) where x ~ r/rs is the radius normalized by a 
scale radius, the best-fit slopes imply that the mass loss 
occurs in the subhalo outskirts beyond the scale radius 
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Fig. 13. — Subhalo mass function spanning two orders of subhalo 
masses. Red solid and blue dotted lines arc the best-fit power-law 
model and Schcchtcr function, respectively. The best-fit powers 
are in remarkable agreement with CDM predictions. Green dashed 
lines are the mass function for spurious peaks. The thick and 
thin dashed lines are the best-fit and the 68% C.L. uncertainty, 
respectively. 

as long as the internal structure dose not change during 
moving in the host halo. 

6.3. Radial Dependence of Subhalo Properties 

Subhalos captured by more massive halos are subject 
to dynamical friction, loses their angular momentum and 
subsequently falls inward the center. Simultaneously, 
their masses arc reduced by the tidal force which in- 
creases with increasing the radius from the cluster center. 
The subhalos in a central region have been affected by 
the tidal field for a longer time than those in outskirts. 
It is thus expected that the subhalo mass and trunca- 
tion radius are an increasing function of cluster radius. 
The survey of subhalos using our wide-field imaging data 
allows us to study the radial dependence of their prop- 
erties. For this purpose, we use the projected position 
of shear-selected subhalos. It is however difficult to con- 
strain a pericenter radius and a line-of-sight position of 
the subhalo. In order to reduce these uncertainties, we 
compute the mean subhalo masses and truncation radii 
derived by the stacked Icnsing analyses for the subsam- 
ple divided by their positions from the cluster center (Sec 
3.3). The mean projected distance of the subsample is es- 
timated as a weighted average of projected distances from 
the cluster center for individual subhalos. The weight 
function is given by the tangential distortion signals with 
respect to the subhalo center. The left and right pan- 
els of Fig. 14 display a clear radial dependence of sub- 
halo masses and truncation radii, as expected from the 
tidal destruction. The subhalo masses and radii grad- 
ually increase out to ~' 1 h^^Mpc and drastically rise 
due to massive subhalo '32' in the outskirts of 1.2 — 
1.6 /i~^Mpc. We fit the subhalo mass profile with a func- 



tional form of log(Afsub/Mpivot) = A + Blog(r/rpivot) 
where Afpivot = 10^^ h~^MQ and rpivot = 1 h^^Mpc and 
obtain, A = 2.66±0.08 and B = 1.45±0.12. The best-fit 
parameters for the inner three bins are A — 1.86 ± 0.23 
and B = 0.55 ± 0.26. The subhalo size profile is fitted 
with log(rt/rt, pivot) = A + Blog{r/rpWot) and rt^pWot = 
1 hr^kpc. The best-fit parameters for all and inner three 
bins are A = 4.95 ± 0.04 and B = 1.18 ± 0.08, and 
A = 4.07 ± 0.26 and B = 0.38 ± 0.24, respectively. The 
best-fit slopes for the mass and truncation radius are sig- 
nificantly changed by the presence of the massive subhalo 
in outskirts (60f <r < 80.'). 

The tidal radius of subhalos is in general defined by 
a competition between the differential tidal forces of 
the host halo potential and the acceleration toward the 
subhalos. Equivalently, this condition can be rewrit- 
ten a balance between the average density of subhalos 



and the host halo, psub 



VPn 



where rj is an ef- 



ficient factor. For instance, 77 = 3 for a point mass 
case on a circular orbit, 77 = 2 for the Roche limit, and 
77 = 2 — dlnMnjain/dlnR for the case of extended mass 
profiles of subhalos and host halo, (e.g. Tormen et al. 
1998; Taylor & Babul 2004; Gan et al. 2010) in a lin- 
ear regime of Mgub/A^fmain ^ 1 and rt/R -^ 1, where 
R is the pericenter radius from the cluster center. A 
minimum subhalo size is determined by a pericenter ra- 
dius. Although we cannot constrain pericenter radius 
of subhalos from the current position of subhalos, it is 
interesting to compare the density ratio, 77, with these 
trial approximation. We calculate the mean densities 
for subhalos, p^nh = Msuh/rt, and for the host halo, 
Pmain = -^%Fw(< R) / R^ ^ whcrc A'/nfw IS a sphcrical 
NFW mass enclosed within the three-dimensional radius. 
Here, we use the best-fit NFW model and the projected 
cluster-centric radius for subhalos. The mean densities 
for subhalos is higher than those for the cluster mass. 
The density ratios, 77 ^^ 10 — 40 for shear-selected sub- 
halos (Sec. 3.3) and ^ 10 for luminous galaxies (Sec. 4) 
are comparable to each other, but higher than expected 
by the linear regime assuming that the current position 
is the pericenter. If the discrepancy could be explained 
by a difference of the positions, shear-selected subhalos 
and luminous galaxies were located inward. On the other 
hand, subhalo mass implies that the detected subhalos 
are remnants of group-scale structure. Large and massive 
subhalos would not be applied by the linear regime. 

6.4. Surface Number Density for subhalos 

It is important to estimate a surface number density 
profile of subhalos for understanding the evolution of 
subhalos in the cluster. Fig. 15 shows the surface num- 
ber density profile of subhalos normalized by the mean 
surface density. The errors for the surface number den- 
sity are assumed to be Poisson noise. It is clear that 
the surface number density increases with decreasing 
the cluster-centric radius. For comparisons, we compute 
the surface number density profiles of member galax- 



ies of which luminosities are L,/ > 10^°/i ^L, 



lO^/i-^L, 



0,i' 



< L,v < 10i°/i-^L, 



0,i' 



and 



©,»'■ 



Both profiles are 
similar to that of subhalos. To quantify these distribu- 
tions, we assume the spherical symmetric NFW distri- 
bution for the subhalos and member galaxies and fit the 
profiles. The surface number density profile is specified 
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Fig. 14. — The mass (left) and truncation radius (right) profiles of subhalos as a function of the projected cluster-centric radius. The 
black circles represent the masses and truncation radii of subhalos for stacked wcak-lcnsing analysis for lensing-selected subhalos, showing 
that they monotonically increase with increasing the radius. 



by three parameters of the normahzation, the concen- 
tration and the virial radius. We here use the virial 
radius determined by the tangential shear fit (Sec 5). 
The best-fit concentrations are Cvir= 5.73 ± 4.46 for 
subhalos, c^-„= 5.97 ± 3.28 for luminous galaxies and 
Cvir= 5.35±1.00 for less luminous galaxies. Interestingly, 
all best-fit concentrations agree with each other. We also 
compute the surface mass density profile for the total 
mass from the best-fit NFW parameters (Sec 5). Here, 
the normalization is set to be the surface density of the 
subhalos at r2oo- The best- fit concentrations for the sub- 
halos and luminous member galaxies do not differ from 
that for the main mass, while the less luminous galaxies 
are more centrally concentrated. Recent numerical simu- 
lations (e.g. Ghigna et al. 2000; Diemand ct al. 2004; De 
Lucia et al. 2004; Gao et al. 2004b, a; Nagai & Kravtsov 
2005; Gao ct al. 2012) and analytical models (e.g. Taylor 
& Babul 2005b; Zentner et al. 2005) have shown that the 
radial distribution of subhalos is less concentrated than 
that of the total mass, because the subhalos lose their 
masses more efficiently in the inner regions of the main 
halo. The local surface density of subhalos is noisy and 
there is an uncertainty of their positions along line-of- 
sight. Coma cluster contains famous NGC4839 group in 
the southwest central region (Figure 3), the presence of 
subhalos in the southwest from the cluster center might 
significantly affect the surface profile. We thus need fur- 
ther studies to assess the subhalo distribution and com- 
pare with dark matter distribution of the main cluster, 
using a large sample of clusters, especially nearby clus- 
ters. 

6.5. Mass-to-light Ratio 

The evolution of galaxies are profoundly affected by 
their surrounding environments, such as dark matter ha- 
los. The environmental processes (e.g. Boselli & Gavazzi 
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Fig. 15. — The surface density (back circles) profile of subhalos. 
The best-fit profile is shown by a blue dotted line. For compari- 
son, a red solid line shows the surface mass density for the total 
mass of the main cluster, normalized with the number density at 
''200 • Red dashed lines arc the 68% confidence uncertainties for 
the total mass. Green diamonds and magenta squares denote the 
surface mass densities for member galaxies of which luminosities 
are L^/ > lO^^h-^Lp^A' and l0^h-^Lfr,.i' < L^, < 10^°h-'^Lr. 
respectively. 



^0,i' 



^Q,i' 



^0,i'; 



2006) in an overdensity region like groups or clusters con- 
sume cold gasses in galaxies through star-bursts trig- 
gered by mergers, tidal interaction with other galax- 
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ies, and ram-prcssurc stripping by the gas, and leads 
to halt the star formation. Based on hierarchical struc- 
ture formation scenario, some cluster galaxies have spent 
long time in group scale environments, before being cap- 
tured in their current host halo. These group scale en- 
vironment would play an important role in the evolu- 
tion for these galaxies, rather than by cluster environ- 
ment (Zabludoff & Mulchaey 1998). On the other hand, 
subhalo masses and sizes depend on their initial proper- 
ties, infall epochs and subsequent evolution in the clus- 
ter halo. Therefore, a study of galaxy-dark matter con- 
nection would provide insights into how galaxies form 
and evolve in different mass properties. We here com- 
pare two independent quantities of weak-lensing masses 
of subhalos and luminosities for associated galaxies. It 
is well established observationally that there is a corre- 
lation between luminosity, velocity dispersion and scale 
length for early type galaxies, so-called the fundamental 
plane (Djorgovski & Davis 1987; Dressier et al. 1987). 
These studies estimate dynamical mass tracing the grav- 
itational potential of a galaxy at small scales. Weak- 
lensing analysis measures the mass of dark matter of 
which distribution is extended beyond galaxies. Thus, 
our approaches are complementary to previous studies 
(e.g. Cappehari et al. 2006; Cody et al. 2009). We com- 
pile two stacked lensing results of shear-selected subha- 
los divided by individual masses (TNFW;Sec. 3.3) and 
luminous member galaxies (Sec. 4). Here, we use the 
truncated mass rather than the NFW model. The lu- 
minosity of each sample is estimated by an average over 
all galaxies associated with subhalos with a weight of 
the tangential distortion signals. The left panel of Fig. 
16, we plot the correlation between subhalo masses and 
those galaxy luminosities. The luminosity ranges be- 
tween 10^° — lO^^/i^^i©, indicating that member galax- 
ies in our subsamples are mainly composed of elliptical 
galaxies. The mass increases with increasing luminosity. 
To quantify this trend, we fit with log(Msub/-Mpivot) = 
A-hBlog(Li//Lpivot), where ipWot = lO^^/i-^L©. The 
best-fit slope, B = 1.49 ± 0.16, gives a positive slope 
with 8cr level. The normalization is A = —0.15 ± 0.19. 
Data points show large scatter. To understand the scat- 
ter, a further careful study using other parameters of 
galaxies is needed to constrain the fundamental plane 
between subhalo masses and galaxy properties. We con- 
vert the scaling relation into the mass-to-light ratio, 

M/L = 86.lt}^;S(LWlO^°/i-'i0)°-^^±°-^^[/iAf0/i.',o]. 
We compare our analysis to the mass-to-light ratio 
determined by dynamical masses. Cappellari et al. 
(2006) estimated the scaling relation of M/L = (2.35 ± 
0.19)(i//10i"i/,o)''-32±0-06 for the SAURON sample. 
Our normalization is higher by one order than dynami- 
cal estimates. It is likely a difference of mass measure- 
ments, because weak-lensing analysis measure mass of 
dark matter halos extend beyond galaxy scales. The sim- 
ilar results have been reported by van Uitert et al. (2011). 
They showed that the weak-lensing masses within r2oo is 
about ten times larger than dynamical masses. 

The right panel of Fig. 16 shows a clear feature of 
decreasing mass-to-light ratio towards the cluster cen- 
ter, similarly to the radial dependence of the mass and 
truncation radius (Fig. 14). Here, we use the TNFW 
mass for subsamples divided by the projected distance 



from the cluster center (Sec. 3.3). The luminosity is es- 
timated by the same way as for subsamples divided by 
mass bins. We fit the form of log(Af/L/(M/L)pivot) = 
A + _B log(r/rpivot) to quantify the radial dependence, 
where (M/i)pivot = IHMq/Lq and rpivot = 1 h'^Mpc. 
We obtain A = 5. 76^0 [l and B = 1.35lJ];^^. The 
upper panel shows the mass-to-light ratio of subhalos 
normalized by the cluster mass-to-light ratio, M/L = 
337.4lg2°5^/i [Mo/Li'^Q], within the virial radius derived 
from the tangential fit. The mass-to-light ratios for sub- 
halos in the outskirts (0.7 ^ 0/Ovir ^ 1) are close to the 
unity, while the ratios in the central region accounts for 
17%. This feature is understood by a scenario that the 
dark matter subhalos arc more subjected to the mass 
loss due to tidal truncation than luminous galaxies which 
tend to be in central region for subhalos. Furthermore, 
since the mean luminosity increases towards the clus- 
ter center, it would be also associated with the galaxy 
evolution. Similar trends are suggested by numerical 
simulations (Springel et al. 2001). They found the the 
median mass-to-light ratio gradually increases out to 
1.5 — 2 /i~^Mpc and is saturated beyond. Gao et al. 
(2004a) also found a similar result that the median mass- 
to- light ratio for subhalos is saturated beyond r2oo- The 
saturated values (Springel et al. 2001) are ~ 15 — 20% of 
the cluster mass-to-light ratio in the B-band. The dis- 
crepancy in the ratios might be caused by a difference 
of subhalos mass range of their simulations and our cat- 
alog, because we used shear-selected subhalos and they 
selected cluster galaxies and associated subhalos. 

This study suggests that the mass-to-light ratio for 
cluster subhalos depends on both the luminosity (mass) 
and the cluster-centric radius. To derive more robust 
conclusions, we need a further study using a large sam- 
ple of nearby cluster lensing analysis. However, our re- 
sult suggests that an assumption of the constant scaling 
relation between the mass and the luminosity gives the 
systematic bias on the mass measurement and their sta- 
tistical properties. 

6.6. Future Studies 

We present the direct observation of the dark mat- 
ter subhalo mass function by weak gravitational lensing 
analysis, which is the first evidence of the consistency 
with CDM predictions on cluster sub-scales. It is thus 
an important step towards studying subhalo mass func- 
tion and properties with a large sample of clusters to 
make stringent tests of the nature of dark matter and 
the details of structure formation. 

Although the subhalo mass function is well described 
by the single power law and the Schechter function, it is 
difficult to discriminate between two functions because 
the abundance of high mass subhalos is low. The sub- 
halo mass function stacked over a large sample of clusters 
will enable us to determine more robustly the functional 
form. Furthermore, the shape of the mass function has 
a characteristic feature depending on the masses of the 
other components of dark matter if any. Thus, the sub- 
halo mass function allows us to constrain the nature of 
dark matter and the structure formation. Hierarchical 
structure formation predicts that subhalo mass function 
depends on host halo mass (e.g. Gao et al. 2004b; van den 
Bosch et al. 2005; Shaw et al. 2006). Less massive halos 
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Fig. 16. — Left panel: The mass and luminosity relation for shear-selected subhalos (black circles) and luminous member galaxies (blue 
squares). The mass of shear-selected subhalos is from stacked lensing results of subsamples divided by model-independent masses. The red 
solid and dashed lines arc the best-fit relation and 68% confidence uncertainties., respectively. Right panel: The radial dependence of the 
mass-to-light ratio obtained from the stacked lensing results for subhalo subsamples divided by their cluster-centric radius. Top panel is 
the mass-to-light ratio normalized by the cluster mass-to-light ratio, as a function of the radius normalized by cluster virial radius (rvir or 
'■200 )■ 



first form at higher redshifts where the mean background 
mass density is higher. Subhalos captured by less massive 
halos are efficiently lost their masses by high density envi- 
ronment. These subhalos are furthermore exposed to de- 
struction for a longer time. Less massive halos arc there- 
fore expected to contain less subhalos. To investigate the 
parent mass dependence, an increase of sample of clus- 
ters is essential. The systematic survey for nearby clus- 
ters for these properties will increase the total number 
of dark matter subhalos to the order of a few hundreds 
or more and improve the statistical accuracy. Further- 
more, finer weak- lensing resolution of nearby clusters will 
enables us to conduct the principal component analyses 
of properties of dark matter halos/subhalos (e.g. Jeeson- 
Danicl et al. 2011; Skibba & Maccio 2011; Wong & Tay- 
lor 2012). Analytical models such as the halo occupation 
distribution (e.g. Scljak 2000; Cooray & Sheth 2002), the 
abundance matching (c.g Vale & Ostriker 2004), and the 
conditional luminosity function (e.g van den Bosch et al. 
2003) would be helpful for understand the galaxy-dark 
matter connection. 

The near future multi-wavelength study of subhalos 
will give us direct and important information on the 
long standing problem of the interplay between dark halo 
and baryons of member galaxies and gasses. We investi- 
gated the correlations between galaxy luminosities, sub- 
halo masses and their projected radius from the cluster 
center. It would be also interesting to investigate the 
correlation with ages, star formation rates and specific 
star formation rates of galaxies. Smith et al. (2012) have 
shown an anisotropic spatial distribution of the galaxies 
age and found that the old population of galaxies dis- 
tribute around the subhalo '17'. It might suggest that 



stellar population properties would vary from subhalo- 
to-subhalo because some cluster galaxies spent long time 
in group scale environments before being captured by the 
cluster. Further systematic studies using other data-sets, 
such as stellar masses, star formation rate, specific star 
formation rate and galaxy types, will provide us with an 
information of cluster galaxy evolution. 

The ROS AT X-ra,y surface brightness distribution (Fig 
5) shows that all of shear-selected subhalos do not con- 
tain X-ray extended structures. It could be understood 
in terms of observational and/or physical effects. First, 
relatively small X-ray sources are unresolved by the poor 
point-spread-function of the ROS AT. Recent observation 
(Andrade-Santos et al. 2013) using Chandra and XMM- 
Newton satellites with high resolutions resolves three X- 
ray subhalos in the central region. They are associated 
with subhalos '21', '23' and '24' , respectively. Subhalo 
masses were estimated under a assumption of hydrostatic 
equilibrium and in good agreement with our mass esti- 
mates. However, X-ray subhalos in other shear-selected 
subhalos are not found even with high resolution data, 
which is due to various physical processes on gas initially 
bound in subhalos. Tormen et al. (2004) has pointed 
out using numerical simulations that the gaseous com- 
ponents, which is collisional matter, are easily destroyed 
by ram pressure stripping and hydrodynamic instabili- 
ties. Accordingly, the lifetime of the X-ray subhalos is 
much shorter than those of dark matter subhalos. Fur- 
thermore, the temperature of the ICM in the central 
region (~ 8keV) is too high to be trapped by subha- 
los. Since X-ray observation of the outskirts requires 
the stable and low X-ray backgrounds, the Suzaku X- 
ray satellite has a great advantage to search gas com- 
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poncnts associated with outskirts subhalos (e.g. Kawa- 
harada et al. 2010; Walker et al. 2012; Sato et al. 2012; 
Ichikawa et al. 2013) rather than the ROSAT. Simionescu 
et al. (2013) has measured the temperature profile out to 
70' ~ 1.4 /i~^Mpc using the Suzaku. The temperature 
in the outskirts drops down to 2keV. The sound veloc- 
ity of this temperature, Cg ^ 720(fcBT/2 keV)^/^ km/s, 
is lower than the escape velocity, Wosc ^ 1500 km/s, ex- 
pected from the most massive subhalo '32'. The enhance 
of gas distribution is thus expected to be detected in this 
region. However, the Suzaku pointings do not fully cover 
the whole area of the cluster. The X-ray follow-up obser- 
vation of the subhalo regions provides us with an impor- 
tant information on the gas evolution and the interplay 
with subhalos (Tozzi & Norman 2001), and resolves a 
possible systematics on the temperature measurement by 
gas clumpy structures in outskirts (Nagai & Lau 2011). 
The thermal Sunyaev Zeldovich effect (SZE) observation 
with different sensitivity from the X-ray is also powerful 
for the gas study. Indeed, Planck Collaboration et al. 
(2012) has shown that the SZE map with FHWM= 10.' 
Planck Collaboration et al. (2012) detected the excess 
flux around NGC 4839 group (subhalo '9'), similar to 
i?05'Ar X-ray image. 

7. CONCLUSIONS 

We conducted weak-lensing survey of subhalos in the 
nearby Coma cluster, with 18 pointing Subaru/Suprime- 
Cam observations, covering 4 square degree and convinc- 
ingly detected 32 subhalos down to the order of 10""^ of 
the virial mass. We quantified systematic issues relevant 
to lensing signals from the large-scale structure behind 
the cluster and a probability of spurious peaks. Our find- 
ings are summarized as follows: 

• Weak-lensing analysis for the nearby Coma cluster 
offers three important advantages to study cluster 
subhalos. First, a large apparent size of subha- 
los enables us to easily resolve less massive subha- 
los. Second, a large apparent area covering sub- 
halos provides a correspondingly large number of 
background galaxies, which leads to low statisti- 
cal errors, compensates low lensing efficiency and 
achieve high signal-to-noise ratio. Third, subha- 
los mass measurements do not suffer from the con- 
taminations in lensing signals from the main clus- 
ter and other subhalos, because subhalos are well 
separated. On the other hand, the analysis has a 
disadvantage to distinguish cluster subhalos from 
background group structures, though we took into 
account the LSS lensing effect in our analysis. A 
spectroscopic or photometric redshifts of galaxies 
is essential to make a secure selection of subhalos. 
A difference of the models of subhalos and back- 
ground groups/clusters is also helpful to assess the 
observed lensing signals. 

• Reconstructed mass maps are associated with the 
projected distributions of member galaxies and LSS 
lensing model at ~ 7 — 14cr, suggesting that ob- 
served shear catalog contains full information of the 
mass structure of and behind the cluster. 

• We discover 32 cluster subhalos by applying thresh- 
olds of peaks appeared in the mass maps. We es- 



timate the model-independent projected masses of 
subhalos, ~ 2— 50x 10^^ H'^Mq, where the smooth 
mass component for the main cluster is subtracted. 

• Staking lensing analysis for samples divided by 
subhalo masses and cluster-centric radii shows a 
sharply truncated profile. The profile is propor- 
tional to (7+ ex 0~^ outside the truncation radii. 
This feature is well described by truncated NEW 
(TNEW) profile rather than the universal NEW 
profile without any truncations, as expected from 
tidal destruction model. As for two subsamples 
of the most and least massive subhalos, the NEW 
model is strongly disfavored. The stacked lensing 
masses are consistent with model-independent ones 
for individual subhalos. 

• The cluster galaxy-galaxy lensing analysis for lu- 
minous member galaxies shows a curved tangential 
shear profile which is well fitted by the NEW model 
or the summation of the truncated NEW model 
considering Gaussian distribution of the truncation 
radius. However, the NEW model is unlikely for 
subhalo models, because the best-fit virial radius 
is too large as the subhalo model. 

• Subhalo mass function dn/d In M^ub is computed 
taking into account systematics of subhalo selec- 
tions. The mass function in the range of two or- 
ders of magnitude in mass is well described by 
single power law ex M~^'^ or the Schechter func- 
tion ex M~J^_^ exp(— A/sub/A^*)- These best-fit slopes 
a = l.lltaag and (3 = 0.94to:29 are in remarkably 
agreements with CDM predictions ~ 0.9 — 1.0 (e.g. 
Giocoh et al. 2010; Gao et al. 2012). This is the 
first evidence of the consistency with CDM predic- 
tions on cluster sub-scales. 

• The subhalo masses, truncation radii and mass-to- 
light ratios decrease toward the cluster center, as 
expected from tidal destruction. The galaxy lumi- 
nosities associated with subhalos depends on both 
their mass and the cluster-centric radius. 

• The tangential distortion signals, g+, in the range 
of ^ 0.02 — 2 /i^^Mpc show a complex structure 
which is well described by three mass components 
of the smooth mass distribution of the NEW model, 
subhalos and LSS lensing model. Although the 
lensing signals are one order lower than those for 
clusters at intermediate redshifts, z ^ 0.2 (Ok- 
abe et al. 2010b), the total signal-to-noise ratio, 
S/N ~ 13.3, is comparable to or higher than be- 
cause of a correspondingly large number of back- 
ground galaxies ('--^ 6 x 10^). The signal-to- noise 
ratios for subhalos and LSS lensing models are 
S/N ~ 4.4 and ~ 1.3, respectively. The 45 de- 
gree rotated component, gx, is consistent with a 
null signal. 
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APPENDIX 

MAP MAKING 

We reconstruct the projected mass distribution following Kaiser & Squires (1993). As described in detail in Okabe 
& Umetsu (2008), we pixelize the reduce shear into a regular grid with a Gaussian smoothing of G{9) oc exp[— 0^/0^]. 

The resolution of maps is defined by FWHM = 2\/\n29g. The smoothed shear pattern at an angular position is 
estimated as 

where Wg and ga^i are the statistical weight of eq. (3) and the reduced shear of the ith galaxy, respectively. The error 
variance for the smoothed shear (Al) is given by 

(E,G(«-«iH,i)' 

Then, the smoothed shear field (Al) is inverted with the kernel (Kaiser & Squires 1993) in Fourier space to obtain the 
projected mass distribution, k{6). Here, we assume the weak-limit of ga = 7q/(1 ~ ^^) ~ la- We also compute the map 
of the significance level for mass reconstruction, v{6) = n/a^, with the mass reconstruction error a^id)- We also make 
maps of luminosity {l{6)) and number density {n(d)) for member galaxies defined in Sec. 2.5 and the convergence field 
(«^LSs(^)) of LSS lensing in Sec. 2.6 with a statistical weight of Wg^i = 1. 

MODEL-INDEPENDENT MASS MEASUREMENT 

A parameter-free estimation of subhalos is given by the aperture-densitomctry, or so-called the Cc-statistics (Clowe 
et al. 2000). The projected mass, Mq^{< 9), is given by 

M^S< e)=7r0'ScrCc(e,ei„n,^out), (Bl) 

Cc{9; 9inn, 9out) = k(< 9) — K{9inn < 9 < 9out) (B2) 

d\n9'{i+{9)) + J I °" d\n9'{^+{9)) 

where ^jnn and 0out are the inner and outer radii of the background annulus. The (7+) is an azimuthal average of the 
tangential component of the gravitational shear, which we take (7+(0)) « {g-t-[9)) in the weak limit. The uncertainty 
in Qc at 9i is estimated as 

where TVinn and A'out are the indices of the discretized radial bins corresponding to the radii of ^jnn and 9out in Eq. (Bl), 
respectively. An error covariance of C,c between each bin is given by (Jij ~ dji = a'^ for 9i < 9j . The signal-to- noise 
ratio of the radial profile, which is complementary information to that of peaks in mass maps, is computed by 

S|N=\Y^M^^,,Vr:HI^^ (B4) 

where K-7 is the inverse of the covariance matrix. 

'■J 

To quantify subhalos' masses, we estimate an average projected mass taking into account the error covariance matrix, 
as follows, 

(B5) 
(B6) 

(B7) 

(B8) 
where TVgi and Ns2 are the indices of the discretized radii where Mr^ profile is saturated. 
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MASS MODELS 

Mass models for cluster halos and sublialos are summarized here. Numerical simulations, based on the CDM model 
of structure formation, predicts that dark matter halos spanning a wide mass range can be described by a universal 
mass density profile (Navarro et al. 1996, 1997). In this paper, we refer to the density profiles as the NFW profile 
which is expressed in the form of 

PNFw(r) = \ — -^-^, (CI) 

{r/rs){l + r/rsy 

where ps is the central density parameter and r,, is the scale radius. The density profile has inner and outer slopes of 
— 1 and —3, respectively. The three-dimensional spherical masses, A/a, enclosed by the radius, ta, inside of which the 
mean density is A times the critical mass density, pcr{z), at the redshift, z, is given by 



Mnfw(< ^a) 



Anps 



ln(l + CA)- Y 



CA 



CA 



(C2) 



The NFW profile is specified by two parameters of Ma and the halo concentration ca — rA/rg. 

Subhalo sizes are determined by the strong tidal field of the main cluster halo. The mass density outside the 
instantaneous tidal radius of subhalos drastically decreases and close to zero due to tidal stripping (e.g. Tormen et al. 
1998; Hayashi ct al. 2003; Oguri & Lee 2004; Taylor & Babul 2004; van den Bosch et al. 2005). Therefore, the mass 
model for subhalos requires an additional parameter of truncation radius, r^. Wc consider two models of the truncation 
model for NFW and SIS models described above. We refer to these truncation models as truncated NFW (TNFW; 
Takada & Jain 2003; Hamana et al. 2004). The interior mass density profile for the TNFW model follows the NFW 
model, but is zero outside the truncation radius. 



Ptnfw (r)^ Pnfw {r <rt), 

= {r>rt). 



(C3) 



This is an extreme case of the truncation model. The TNFW model (Okabe et al. 2010a) is specified by three 
parameters of subhalo mass (Mt), a truncation radius (r^) and a concentration (ct). A slope of lensing profile for 
the TNFW model drastically changes outside the truncation radius and behaves as a point source (ex 9^^). In model 
fitting, the truncation radius and the concentration are sensitive to this break and the inner profile, respectively. The 
subhalo mass is also determined by the distortion signal at the truncation radius. Since the mass densities outside the 
truncation radius is zero, the three- and two- dimensional masses within the truncation radius yield the exactly same 
value (MsD = M2b)- 

When we stack tangential distortion profiles for subhalos with different properties, the break in the distortion profile 
is smooth due to their intrinsic distribution. Addition to the TNFW model, we compute the model taking into 
account a distribution function of the truncation radius (TNFWProb). The probability function is assumed to be 
Gaussian, p{rt) = exp(— (r^ — (rt))^/2CT^J/-y27rCT^, where (r^) and a^ are the average and the standard error for the 

truncation radius, respectively. The mean lensing signal is expressed in terms of j p(j't)g\. drt- We don't assume 

a distribution on subhalo mass but instead estimate the mean subhalo mass, (Mt), which is sensitive to a lensing signal 

at (rt). 
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